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Abstract. The fact that one can associate thermodynamic properties with horizons 
brings together principles of quantum theory, gravitation and thermodynamics and 
possibly offers a window to the nature of quantum geometry. This review discusses 
certain aspects of this topic concentrating on new insights gained from some recent 
work. After a brief introduction of the overall perspective. Sections 2 and 3 provide 
the pedagogical background on the geometrical features of bifurcation horizons, 
path integral derivation of horizon temperature, black hole evaporation, structure 
of Lanczos-Lovelock models, the concept of Noether charge and its relation to 
horizon entropy. Section 4 discusses several conceptual issues introduced by the 
existence of temperature and entropy of the horizons. In Section 5 we take up the 
connection between horizon thermodynamics and gravitational dynamics and describe 
several peculiar features which have no simple interpretation in the conventional 
approach. The next two sections describe the recent progress achieved in an alternative 
perspective of gravity. In Section 6 we provide a thermodynamic interpretation of the 
field equations of gravity in any diffeomorphism invariant theory and in Section 7 we 
obtain the field equations of gravity from an entropy maximization principle. The last 
section provides a summary. 
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1. Introduction and perspective 

Soon after Einstein developed the gravitational field equations, Schwarzschild found the 
simplest exact solution to these equations, describing a spherically symmetric spacetime. 
When expressed in a natural coordinate system which makes the symmetries of the 
solution obvious, it leads to the line interval: 



where /(r) = 1 — (r^/r) with = 2GM/c^ = 2M, in units with G = c = 1. It 
was immediately noticed that this metric exhibits a curious pathology. One of the 
metric coefficients, gu vanished on a surface "H, of finite area 47rr^, given by r = r^, 
while another metric coefficient grr diverged on the same surface. After some initial 
confusion, it was realized that the singular behaviour of the metric is due to bad choice 
of coordinates and that the spacetime geometry is well behaved at r = r^. However, 
the surface H acts as a horizon blocking the propagation of information from the region 
r < to the region r > Vg. This leads to several new features in the theory, many 
of which, even after decades of investigation, defies a complete understanding. The 
most important amongst them is the relationship between physics involving horizons 
and thermodynamics. 

This connection, between black hole dynamics involving the horizon and the laws of 
thermodynamics, became apparent as a result of the research in early 1970s. Hawking 
proved [1] that in any classical process involving the black holes, the sum of the areas 
of a black hole horizons cannot decrease which, with hindsight, is reminiscent of the 
behaviour of entropy in classical thermodynamics. This connection was exploited by 
Bekenstein in his response to a conundrum raised by John Wheeler. Wheeler pointed out 
that an external observer can drop material with non-zero entropy into the inaccessible 
region beyond the horizon thereby reducing the entropy accessible to outside observers. J 
Faced with this difficulty, Bekenstein came up with the idea that the black hole horizon 
should be attributed an entropy which is proportional to its area [3, 4, 5]. It was also 
realized around this time that one can formulate four laws of black hole dynamics in a 
manner analogous to the laws of thermodynamics [6] . In particular, if a physical process 
(say, dropping of small amount of matter into the Schwarzschild black hole) changes the 
mass of a black hole by 5M and the area of the event horizon by 5A, then it can be 
proved that 



where k = M/(2M)^ = 1/4M is called the surface gravity of the horizon. This suggests 
an analogy with the thermodynamic law 6E = T6S with S oc A and T oc k. However, 

I John Wheeler posed this as a question to Bekenstein: What happens if you mix cold and hot tea 
and pour it down a horizon, erasing all traces of "crime" in increasing the entropy of the world? This 
is based on what Wheeler told me in 1985, from his recollection of events; it is also mentioned in his 
book, see page 221 of Ref. [2]. I have heard somewhat different versions from other sources. 



ds^ = -f{r)dt^ + 



dr"^ 



+ r2(c/^2 + sin2 Odcj)'^) 



(1) 




(2) 
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classical considerations alone cannot determine the proportionality constants. (We will 
see later that quatum mechanical considerations suggest S = {A/ 4) and T = h/2ti, 
which is indicated in the second equality in Eq. (2).) 

In spite of this, Bekenstein's idea did not find favour with the community 
immediately; the fact that laws of black hole dynamics have an uncanny similarity 
with the laws of thermodynamics was initially considered to be only a curiosity. (For a 
taste of history, see e.g [7].) The key objection at that time was the following. If black 
holes possess entropy as well as energy (which they do), then they must have a non-zero 
temperature and must radiate — which seemed to contradict the view that nothing can 
escape the black hole horizon. Investigations by Hawking, however, led to the discovery 
that a non-zero temperature should be attributed to the black hole horizon [8] . He found 
that black holes, formed by the collapse of matter, will radiate particles with a thermal 
spectrum at late times, as detected by a stationary observer at large distances. This 
result, obtained from the study of quantum field theory in the black hole spacetime, 
showed that one can consistently attribute to the black hole horizon an entropy and 
temperature. 

An immediate question that arises is whether this entropy is the same as the "usual 
entropy". If so, one should be able to show that, for any processes involving matter 
and black holes, we must have d{SBH + Smauer)/dt > which goes under the name 
generalized second law (GSL). One simple example in which the area (and thus the 
entropy) of the black hole decreases is in the emission of Hawking radiation itself; but 
the GSL holds since the thermal radiation produced in the process has entropy. It is 
generally believed that GSL always holds though a completely general proof is difficult 
to obtain. Several thought experiments, when analyzed properly, uphold this law (see, 
for example, Ref. [9]) and a proof is possible under different sets of assumptions [10]. 
All these suggest that the area of the black hole corresponds to an entropy which is 
same as the "usual entropy" . 

These ideas can be extended to black hole solutions in more general theories than 
just Einstein's gravity. The temperature T can be determined by techniques like, for 
example, analytic continuation to imaginary time (see Sec. 3.2) which depends only 
on the metric and not on the field equations which led to the metric. But the concept 
of entropy needs to generalized in these models and will no longer be one quarter of 
the area of horizon. This could be done by using the first law of black hole dynamics 
itself, say, in the form TdS = dM. Since the temperature is known, this equation 
can be integrated to determine 5*. This was done by Wald and it turns out that the 
entropy can be related to a conserved charge called Noether charge which arises from 
the diffeomorphism invariance of the theory [11]. Thus, the notions of entropy and 
temperature can be attributed to black hole solutions in a wide class of theories. 

This raises the question: What are the degrees of freedom responsible for the black 
hole entropy? There have been several attempts in the literature to answer this question 
both with and without inputs from quantum gravity models. A statistical mechanics 
derivation of entropy was originally attempted in [12]; the entropy has been interpreted 
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as the logarithm of: (a) the number of ways in which black hole might have been formed 
[4, 13]; (b) the number of internal black hole states consistent with a single black hole 
exterior [14] and (c) the number of horizon quantum states [2, 15, 16]. There are also 
other ideas which are more formal and geometrical [17, 18, 19, 20], or based on thermo- 
field theory [21, 22, 23] just to name two possibilities. In addition, considerable amount 
of work has been done in calculating black hole entropy based on different candidate 
models for quantum gravity. (We will briefly summarize some aspects of these in Sec. 
4.2; more extensive discussion as well as references to original literature can be found 
in the reviews [24, 25]). 

Clearly, there is no agreement in literature as regards the degrees of freedom which 
contribute to black hole entropy and the attempts mentioned above sample the divergent 
views. In fact, once the answer is known, it seems fairly easy to come up with very 
imaginative derivations of the result. 

A crucial new dimension is added to this problem when we study horizons which 
are not associated with black holes. Soon after Hawking's discovery of a temperature 
associated with the black hole horizon, it was realized that this result was not confined 
to black holes alone. The study of quantum field theory in any spacetime with a horizon 
showed that all horizons possess temperatures [26, 27, 28]. In particular, an observer 
who is accelerating through the vacuum state in fiat spacetime perceives a horizon and 
will attribute to it [29] a temperature T = k/27i proportional to her acceleration k. 
(For a review, see Refs.[30, 31, 32, 33, 34, 35, 36].) The situation regarding entropy 
— especially whether one should attribute entropy to all horizons — remains unclear. 
In fact, many of the attempts to interpret black hole entropy mentioned earlier cannot 
be generalized to interpret the entropy of other horizons. So a unified understanding 
of horizon thermodynamics encompassing both entropy and temperature still remains 
elusive — which will be one of the key issues we will discuss in detail in this review. 

This is also closely related to the question of whether gravitational dynamics — 
in particular, the field equations — have any relation to the horizon thermodynamics. 
Given a spacetime metric with a horizon (which may or may not be a solution to 
Einstein's equations) one can study quantum field theory in that spacetime and discover 
that the horizon behaves like a black body with a given temperature. At no stage in such 
an analysis do we need to invoke the gravitational field equations. So it is reasonable to 
doubt whether the dynamics of gravity has anything to do with horizon thermodynamics 
and one may — at first — think that there is no connection between the two. 

Several recent investigations have shown, however, that there is indeed a deeper 
connection between gravitational dynamics and horizon thermodynamics (for a recent 
review, see Ref. [37]). For example, studies have shown that: 

• Gravitational field equations in a wide variety of theories, when evaluated on a 
horizon, reduce to a thermodynamic identity TdS = dE + PdV . This result, 
first pointed out in Ref. [38], has now been demonstrated in several cases like the 
stationary axisymmetric horizons and evolving spherically symmetric horizons in 
Einstein gravity, static spherically symmetric horizons and dynamical apparent 
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horizons in Lovelock gravity, three dimensional BTZ black hole horizons, FRW 
cosmological models in various gravity theories and even in the case Horava-Lifshitz 
gravity (see Sec. 5.1 for detailed references). If horizon thermodynamics has no deep 
connection with gravitational dynamics, it is not possible to understand why the 
field equations should encode information about horizon thermodynamics. 

• Gravitational action functionals in a wide class of theories have a a surface term and 
a bulk term. In the conventional approach, we ignore the surface term completely 
(or cancel it with a counter-term) and obtain the field equation from the bulk 
term in the action. Therefore, any solution to the field equation obtained by this 
procedure is logically independent of the nature of the surface term. But when the 
surface term (which was ignored) is evaluated at the horizon that arises in any given 
solution, it gives the entropy of the horizon! Again, this result extends far beyond 
Einstein's theory to situations in which the entropy is not proportional to horizon 
area. This is possible only because there is a specific holographic relationship 
[39, 40, 41] between the surface term and the bulk term which, however, is an 
unexplained feature in the conventional approach to gravitational dynamics. Since 
the surface term has the thermodynamic interpretation as the entropy of horizons, 
and is related holographically to the bulk term, we are again led to suspect an 
indirect connection between spacetime dynamics and horizon thermodynamics. 

Based on these features — which have no explanation in the conventional approach 
— one can argue that there is a conceptual reason to revise our perspective towards 
spacetime (Sec. 6 and 7) and relate horizon thermodynamics with gravitational 
dynamics. This approach should work for a wide class of theories far more general 
than just Einstein gravity. This will be the new insight which we will focus on in this 
review. 

To set the stage for this future discussion, let us briefly describe this approach and 
summarize the conclusions. We begin by examining more closely the implications of the 
existence of temperature for horizons. 

In the study of normal macroscopic systems — like, for example, a solid or a gas — 
one can deduce the existence of microstructure just from the fact that the object can be 
heated. The supply of energy in the form of heat needs to be stored in some form in the 
material which is not possible unless the material has microscopic degrees of freedom. 
This was the insight of Boltzmann which led him to suggest that heat is essentially 
a form of motion of the microscopic constituents of matter. That is, the existence of 
temperature is sufficient for us to infer the existence of microstructure without any 
direct experimental evidence. 

The thermodynamics of the horizon shows that we can actually heat up a spacetime, 
just as one can heat up a solid or a gas. An unorthodox way of doing this would be 
to take some amount of matter and arrange it to collapse and form a black hole. The 
Hawking radiation emitted by the black hole can be used to heat up, say, a pan of 
water just as though the pan was kept inside a microwave oven. In fact the same result 



Thermodynamical Aspects of Gravity: New insights 



6 



can be achieved by just accelerating through the inertial vacuum carrying the pan of 
water which will eventually be heated to a temperature proportional to the acceleration. 
These processes show that the temperatures of the horizons are as "real" as any other 
temperature. Since they arise in a class of hot spacetimes, it follows d la Boltzmann 
that the spacetimes should possess microstructure. 

In the case of a solid or gas, we know the nature of this microstructure from atomic 
and molecular physics. Hence, in principle, we can work out the thermodynamics of 
these systems from the underlying statistical mechanics. This is not possible in the case 
of spacetime because we have no clue about its microstructure. However, one of the 
remarkable features of thermodynamics — in contrast to statistical mechanics — is that 
the thermodynamic description is fairly insensitive to the details of the microstructure 
and can be developed as a fairly broad frame work. For example, a thermodynamic 
identity like TdS = dE + PdV has a universal validity and the information about 
a given system is only encoded in the form of the entropy functional S{E^V). In 
the case of normal materials, this entropy arises because of our coarse graining over 
microscopic degrees of freedom which are not tracked in the dynamical evolution. In 
the case of spacetime, the existence of horizons for a particular class of observers makes it 
mandatory that these observers integrate out degrees of freedom hidden by the horizon. 

To make this notion clearer, let us start from the principle of equivalence which 
allows us to construct local inertial frames (LIF), around any event in an arbitrary 
curved spacetime. Given the LIF, we can next construct a local Rindler frame (LRF) 
by boosting along one of the directions with an acceleration k. The observers at rest in 
the LRF will perceive a patch of null surface in LIF as a horizon T-L with temperature 
K/27r. These local Rindler observers and the freely falling inertial observers will attribute 
different thermodynamical properties to matter in the spacetime. For example, they will 
attribute different temperatures and entropies to the vacuum state as well as excited 
states of matter fields. When some matter with energy 5E moves close to the horizon 
— say, within a few Planck lengths because, formally, it takes infinite Rindler time 
for matter to actually cross l-L — the local Rindler observer will consider it to have 
transfered an entropy 5S = {2tt / k)5E to the horizon degrees of freedom. We will show 
(in Sec. 6) that, when the metric satisfies the field equations of any diffeomorphism 
invariant theory, this transfer of entropy can be given [42] a geometrical interpretation 
as the change in the entropy of the horizon. 

This result allows us to associate an entropy functional with the null surfaces which 
the local Rindler observers perceive as horizons. We can now demand that the sum of 
the horizon entropy and the entropy of matter that flows across the horizons (both as 
perceived by the local Rindler observers), should be an extremum for all observers in 
the spacetime. This leads [43] to a constraint on the geometry of spacetime which can 
be stated, in D = 4, as 

{Gab - STTTatKn" = (3) 
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for all null vectors in the spacetime. The general solution to this equation is given 
by Gab = SvrTab + poQab where po has to be a constant because of the conditions 
VaG""^ = = VaT""^. Hence the thermodynamic principle leads uniquely to Einstein's 
equation with a cosmological constant in 4-dimensions. Notice, however, that Eq. (3) 
has a new symmetry and is invariant [44, 45] under the transformation Tab Tab + ^gab 
which the standard Einstein's theory does not posses. (This has important implications 
for the cosmological constant problem [46] which we will discuss in Sec. 7.5.) In D > 4, 
the same entropy maximization leads to a more general class of theories called Lanczos- 
Lovelock models (see Sec. 3.6). 

We can now remedy another conceptual shortcoming of the conventional approach. 
An unsatisfactory feature of all theories of gravity is that the field equations do not have 
any direct physical interpretation. The lack of an elegant principle which can lead to 
the dynamics of gravity ( "how matter tells spacetime to curve" ) is quite striking when 
we compare this situation with the kinematics of gravity ("how spacetime makes the 
matter move"). The latter can be determined through the principle of equivalence by 
demanding that all freely falling observers, at all events in spacetime, must find that 
the equations of motion for matter reduce to their special relativistic form. 

In the alternative perspective, Eq. (3) arises from our demand that the 
thermodynamic extremum principle should hold for all local Rindler observers. This 
is identical to the manner in which freely falling observers are used to determine how 
gravitational field influences matter. Demanding the validity of special relativistic laws 
for the matter variables, as determined by all the freely falling observers, allows us to 
determine the influence of gravity on matter. Similarly, demanding the maximization 
of entropy of horizons (plus matter), as measured by all local Rindler observers, leads 
to the dynamical equations of gravity. 

In this review, we will examine several aspects of these features and will try to 
provide, in the latter part of the review, a synthesis of ideas which offers an interesting 
new perspective on the nature of gravity that makes the connection between horizon 
thermodynamics and gravitational dynamics obvious and natural. In fact, we will 
show that the thermodynamic underpinning goes far beyond Einstein's theory and 
encompasses a wide class of gravitational theories. 

The review is organized as follows. In Section 2 we will review several features 
of horizons which arise in different contexts in gravitational theories. In particular, 
we will describe some generic features of the spacetimes with horizons which will be 
important in the later discussions. In Section 3 we shall provide a simple derivation 
of the temperature of a (generic) horizon using path integral methods. Section 3.1 
introduces the basic concepts related to path integrals and Section 3.2 applies them to 
a spacetime with horizon to obtain the temperature. Some alternate ways of obtaining 
the temperature of horizons are discussed in Section 3.3. The origin of Hawking 
radiation from matter that collapses to form a black hole is discussed in Section 3.4. 
The next two subsections describe the generalization of these ideas to theories other 
than Einstein's general relativity. The relationship between horizon entropy and the 
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Noether charge is introduced in Section 3.5 and the structure of Lanczos-Lovelock 
models is summarized in Section 3.6. Section 4 discusses several conceptual issues 
raised by the existence of temperature and entropy of the horizons, concentrating on 
the nature of degrees of freedom which contribute to the entropy and the observer 
dependence of the concept of entropy. This discussion is continued in the next section 
where we take up the connection between horizon thermodynamics and gravitational 
dynamics and describe several peculiar features which have no simple interpretation in 
the conventional approach. In Section 6 we provide a thermodynamic interpretation of 
the field equations of gravity in any diffeomorphism invariant theory. This forms the 
basis for the alternative perspective of gravity described in Section 7 in which we obtain 
the field equations of gravity from an entropy maximization principle. The last section 
provides a summary. 

We will use the signature ( — h ++) and units with G = h = c = 1. The Greek 
superscripts and subscripts will run over the spatial coordinates while the Latin letters 
will cover time coordinate as well as spatial coordinates. 



2. Gravity and its horizons 

2.1. The Rindler horizon in flat spacetime 

The simplest context in which a horizon arises for a class of observers occurs in the flat 
spacetime itself. Consider the standard fiat spacetime metric with Cartesian coordinates 
in the X — T plane given by 

ds^ = -dT^ + dX^ + dLl (4) 

where dL^ is the line element in the transverse space. The lines X = ±T divide the 
X — T plane into four quadrants (see Fig. 1) marked the right (7^) and left (£) wedges as 
well as the past {V) and future (J-") of the origin. We now introduce two new coordinates 
[t, I) in place of (T, X) in all the four quadrants through the transformations: 

kT = V2kI sinh(Kt); nX = ±\/2kI cosh(/tt) (5) 

for |X| > \T\ with the positive sign in TZ and negative sign in C and 

kT = ±y/—2hil cosh(Kt); nX = y/—2hil sinh(Kt) (6) 

for |X| < |T| with the positive sign in J-" and negative sign in V. Clearly, / < is used 
in J-" and V. With these transformations, the metric in all the four quadrants can be 
expressed in the form 

ds^ = -2kI de + — + dLl (7) 
2kI 

Figure 1 shows the geometrical features of the coordinate systems from which we see 
that: (a) The coordinate t is timelike and / is spacelike in Eq. (7) only in TZ and C 
where / > with their roles reversed in J-" and V with / < 0. (b) A given value of (t, /) 
corresponds to a pair of points in TZ and C for / > and to a pair of points in J-" and 
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V for / < 0. (c) The surface / = acts as a horizon for observers in TZ. In particular, 
observers who are stationary in the new coordinates with / = constant, x_|_ = constant 
will follow a trajectory X"^ — = 2//k in the X — T plane. These are trajectories of 
observers moving with constant proper acceleration in the inertial frame who perceive a 
horizon at Z = 0. Such observers are usually called Rindler observers and the metric in 
Eq. (7) is called the Rindler metric. (The label N in Fig. 1 corresponds to iV = -\/2k]7[.) 



T 

A 




Figure 1. The global manifold with different coordinate systems in the four quadrants. 
See text for discussion. 



2.2. The Rindler frame as the near-horizon limit 

The Rindler frame will play a crucial role in our future discussions for two reasons. First, 
one can introduce Rindler observers even in curved spacetime in any local region. To 
do this, we first transform to the locally inertial frame (with coordinates T, X) around 
that event and then introduce the local Rindler frame with coordinates (t, /) by the 
transformations in Eq. (5) and Eq. (6). Such a local notion is approximate but will 
prove to be valuable in our future discussions because it can be introduced around any 
event in any curved spacetime. Second, the Rindler (like) transformations work for a 
wide variety of spherically symmetric solutions to gravitational field equations. This 
general class of spacetimes can be expressed by a metric of the form 

ds' = -f{r)de + ^,+dLl (8) 

where the function /(r) has a simple zero at some point r = a with a non-zero first 
derivative /'(a) = 2k. A Taylor series expansion of / near r = a gives / ~ 2k,1 with 
I = r — a. It is, therefore, obvious that near the horizon, located at r = a, all these 
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metrics can be approximated by the Rindler metric in Eq. (7). Hence Rindler metric 
is useful in the study of spacetime near the horizon in several exact solutions. (We are 
assuming that /'(a) 7^ 0; there are certain solutions — called extremal horizons — in 
which this condition is violated and f'{a) = 0; we will not discuss them in this review.) 

In the above analysis we started from a flat spacetime expressed in standard 
inertial coordinates and then introduced the transformation to Rindler coordinates. 
This transformation, in turn, brought in a pathological behaviour for the metric at 
/ = 0. Alternatively, if we were given the metric in Rindler coordinates in the form of 
Eq. (7) we could have used the transformation in Eq. (5) and Eq. (6) to remove the 
pathological behaviour of the metric. In such a process we would have also discovered 
that a given value of {t, I) actually corresponds to a pair of events in the full spacetime 
thereby 'doubling up' the manifold. This process is called analytic extension. 

In the case of metrics given by Eq. (8) the pathology at / = is similar to the 
pathology of the Rindler metric at / = 0. Just as one can eliminate the latter by analytic 
extension, one can also eliminate the singularity at r = a in the metric in Eq. (8) by 
suitable coordinate transformation. Consider for example, the transformations from 
(t, r) to (T, X) by the equations 

/dv 
-TT^ (9) 

This leads to a metric of the form 

ds' = ^2^J_T2^ i-dT^ + dX^) + dLl (10) 

where / needs to be expressed in terms of (T, X) using the coordinate transformations. 
The horizon r = a now gets mapped to X"^ = T^; but it can be shown that the factor 
f /{X"^ — T^) remains finite at the horizon. 

The similarity between the coordinate transformations in Eq. (9) and Eq. (5) is 
obvious. (As in the case of Eq. (5), one can introduce another set of transformations 
to cover the remaining half of the manifold by interchanging sinh and cosh factors.) 
The curves of constant r in the original spherically symmetric metric in Eq. (8) become 
hyperbolas in the T — X plane, just as in the case of transformation from Rindler to 
inertial coordinates. 



2.3. Horizons in static spacetimes 

In the Rindler frame (as well as near the horizon in a curved spacetime), one can 
introduce another coordinate system which often turns out to be useful. This is done 
by transforming from {t,l) to {t,x) where / = {1/2)k,x'^. Then the Rindler metric in 
Eq. (7) reduces to the form 

ds^ = -K^x^dt^ + dx^ + dL\ (11) 

and the coordinate transformation transformation corresponding to Eq. (5) becomes 

T = xsmh^Kt); X = cosh(fi:t) (12) 
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The form of the metric in Eq. (11) also arises in a wide class of static (and stationary, 
though we will not discuss this case) spacetimes with the following properties: (i) 
The metric is static in the given coordinate system, goa = 0,gab(t,^) = gab{'^)', (ii) 
(yfoo(x) = — A^2(x) vanishes on some 2-surface "H defined by the equation A^^ = 0, (iii) 
daN is finite and non zero on T-L and (iv) all other metric components and curvature 
remain finite and regular on Tt. The line element will now be: 

ds^ = -N\x'')dt^ + ^o.i3ix")dx"dx^ (13) 

The comoving observers in this frame have trajectories x = constant, four-velocity 
Ua = —N6^ and four acceleration a* = u^VjU^ = (0, a) which has the purely spatial 
components = {daN)/N. The unit normal Ua to the N = constant surface is given 
by ria = daN [g'^" d^N dyN)'^^"^ = aa{ai3a^)~^/'^ . A simple computation now shows that 
the normal component of the acceleration a*nj = a'^Ua, 'redshifted' by a factor N, has 
the value 

Ar(n„a°) = {g^'^d^NdpNy/'^ = iVa(x) (14) 

where the last equation defines the function a. From our assumptions, it follows that on 
the horizon N = 0, this quantity has a finite limit Na — )■ k; the k is called the surface 
gravity of the horizon. 

These static spacetimes, however, have a more natural coordinate system defined 
in terms of the level surfaces of A^. That is, we transform from the original space 
coordinates x^ in Eq.(13) to the set (A^, y^), A = 2, 3 by treating A^ as one of the spatial 
coordinates — which is always possible locally. The y"^ denotes the two transverse 
coordinates on the A^ = constant surface. The line element in the new coordinates will 
be: 

'''' = * + JN^ + '^-(''f - Iv^X^f - <15) 

where a"* etc. are the components of the acceleration in the new coordinates. The 
original 7 degrees of freedom in (A^, 7^,^) are now reduced to 6 degrees of freedom in 
{a,a^,aAB), because of our choice for goQ. In Eq.(15) the spacetime is described in 
terms of the magnitude of acceleration a, the transverse components and the metric 
a^B on the two surface and maintains the t— independence. The A^ is now merely a 
coordinate and the spacetime geometry is described in terms of (a, a"^, <Jab) all of which 
are, in general, functions of {N,y^). In spherically symmetric spacetimes with horizon, 
for example, we will have a = a{N), a"^ = if we choose y^ = {9, (p). Important features 
of dynamics are usually encoded in the function a{N,y^). Near the A^ — )■ surface, 
Na — K, the surface gravity, and the metric reduces to 

ds^ = -N'dt' + 4^ + dLl ^ -N^de + ^ + dLl (16) 

where the second equality is applicable close to %. This is the same metric as in Eq.(ll) 
if we set N = kx. Therefore a wide class of metrics with horizon can be mapped to the 
Rindler form near the horizon. 
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The form of the metric in Eq. (11) is particularly useful to study the analytic 
continuation to imaginary values of time coordinate. If we denote Te = iT,tE = it, 
then, in the right wedge TZ, the transformations become 

Te = X sin^ntE)', X = xcos{KtE) (17) 

which are just the coordinate transformation from the Cartesian coordinates {Te,X) 
to the polar coordinates {6 = KtE,x) in a two dimensional plane. To avoid a conical 
singularity at the origin, it is necessary that 6 is periodic with period 27r, which — 
in turn — requires to be periodic with period 2tt/k,. We will see later that such a 
periodicity in the imaginary time signals the existence of non-zero temperature. 



2.4- Exponential redshift and thermal power spectrum 

Another generic feature of the horizons we have defined is that they act as surfaces of 
infinite redshift. To see this, consider the redshift of a photon emitted at {t^, Ne,y'^), 
where A'^e is close to the horizon surface "H, and is observed at (t, N, y^). The frequencies 
at emission u(te) and detection u(t) are related by [oj(t)/u{te)] = [Nf,/N]. The radial 
trajectory of the out-going photon is given by ds^ = which integrates to 

dN 1 

t — te = / -T-jY- — In A'e + constant (18) 

where we have approximated the integral by the dominant contribution near N^. = 
This gives oc exp(— /tt), leading to the exponentially redshifted frequency 

(jj{t) oc A'e oc exp(— Kt) (19) 

as detected by an observer at a fixed as a function of t. 

Such an exponential redshift is also closely associated with the emergence of a 
temperature in the presence of a horizon. To see this, let us consider how an observer in 
Rindler frame (or, more generally, in the spherically symmetric frame with the metric 
given by Eq. (8)) will view a monochromatic plane wave moving along the X— axis in 
the inertial frame (or, more generally, in the analytically extended coordinates). Such 
a scalar wave can be represented by (j){T,X) = exp[—iQ{T — X)] with > 0. Any 
other observer who is inertial with respect to the X = constant observer will see this 
as a monochromatic wave, though with a different (Doppler-shifted) frequency. But an 
accelerated observer, at A^ = A^g = constant using her proper time r = Aq^ will see the 
same mode as varying as 

= 0(T(t),X(t)) = exp[inqe-''^] = exp[inqexp -{k/No)t] (20) 

where we have used Eq. (9) and defined q = K~^exp{K^). This is clearly not 
monochromatic and has a frequency which is being exponentially redshifted in time. 
The power spectrum of this wave is given by P{iy) = |/(^)P where /(z/) is the Fourier 
transform of 4>{t) with respect to r: 

/OO 7 /* OO 7 

^ ^/('^)^"^^^ £ [M'^)^-"'^ + B{y)e-] (21) 
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with = /(i/) and -B(i^) = f{—v). Because of the exponential redshift, this power 
spectrum will not vanish for z/ < leading to i? 7^ 0. Evaluating this Fourier transform 
(by changing to the variable r2gexp[— (fi:/Ai'o)r] = z and analytically continuing to Im 
z) one gets: 

This leads to the the remarkable result that the power, per logarithmic band in 
frequency, at negative frequencies is a Planckian at temperature T = {k,/27iNq): 

Hi?(.)P = H/(-)P = ^; /3 = ^ (23) 

Though /(i/) in Eq. (22) depends on Q, the power spectrum |/(i^)P is independent of 
Q; monochromatic plane waves of any frequency (as measured by the freely falling 
observers with X = constant) will appear to have Planckian power spectrum in 
terms of the (negative) frequency u, defined with respect to the proper time of the 
accelerated observer located at = A'o = constant. The scaling of the temperature 
(3~^ oc A^'q'^ oc Is'ooI "'^^^ is precisely what is expected in general relativity for temperature. 
(Similar results also arise in the case of real wave with oc cosr2(T— X); see ref. [47, 48].) 

We saw earlier (see Eq. (18)) that waves propagating from a region near the horizon 
will undergo exponential redshift. An observer detecting this exponentially redshifted 
radiation at late times (t 00), originating from a region close to "H will attribute to 
this radiation a Planckian power spectrum given by Eq. (23). This result lies at the 
foundation of associating a temperature with a horizon. 

The Planck spectrum in Eq. (23) is in terms of the frequency and /3 oc c/k has the 
(correct) dimension of time; no h appears in the result. If we now switch the variable 
to energy and write fSu = {P/h){hv) = {P/h)E, then one can identify a temperature 
ksT = {i%h/2Tic) which scales with h. This "quantum mechanical" origin of temperature 
is superficial because it arises merely because of a change of units from v to E. An 
astronomer measuring frequency rather than photon energy will see the spectrum in 
Eq. (23) as Planckian without any quantum mechanical input. The real role of quantum 
theory is not in the conversion of frequency to energy but in providing the complex wave 
in the inertial frame. It represents the vacuum fluctuations of the quantum field. 



2.5. Field theory near the horizon: Dimensional reduction 



The fact that —Qqq = N"^ — on the horizon leads to several interesting conclusions 
regarding the behaviour of any classical (or quantum) field near the horizon. Consider, 
for example, an interacting scalar field in a background spacetime described by the 
metric in Eq.(15), with the action: 



A 



d'^Xy/^ -da 



dtdNd^y 



+ V 



(24) 



ON 



V 
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where {d±(j))^ denotes the contribution from the derivatives in the transverse directions 
including cross terms of the type {di\f(f)d±(f)). Near = 0, with Na — )■ k, the action 
reduces to the form 




where we have changed variable to ^ defined in Eq. (9) (which behaves as ^ ~ (1/k) In 
near ^ 0) and ignored terms that vanish as — 0. Remarkably enough this action 
represents a two dimensional free field theory in the (t,0 coordinates which has the 
enhanced symmetry of invariance under the conformal transformations Qab P{^i Odab 
(see e.g., Section 3 of [49], [50]). The solutions to the field equations near "H are plane 
waves in the (t,0 coordinates: 

0± = exp[-iu{t ± 0] = A^±^'^/«e"^'^* (26) 

These modes are the same as = expiA where A is the solution to the classical 
Hamilton- Jacobi equation; this equality arises because the divergence of (1/A^) factor 
near the horizon makes the WKB approximation almost exact near the horizon. The 
mathematics involved in this phenomenon is fundamentally the same as the one which 
leads to the "no-hair-theorems" (see, eg., [51]) for the black hole. These solutions possess 
several other symmetry properties which are worth mentioning: 

To begin with, the metric and the solution near "H are invariant under the rescaling 
A^ — )■ AA^, in the sense that this transformation merely adds a phase to 0. This scale 
invariance can also be demonstrated by studying the spatial part of the wave equation 
[52] near Ti, where the equation reduces to a Schrodinger equation for the zero energy 
eigenstate in the potential V^(A^) = —u^/N"^ . This Schrodinger equation has the natural 
scale invariance with respect to N XN which is reflected in our problem. 

Second, the relevant metric ds"^ = —N^dt"^ + {dN / kY in the t — N plane is also 
invariant, up to a conformal factor, to the metric obtained by A^ — )■ p = 1/A^: 

ds^ = -NHe + ^ = -A-p'de + ^) (27) 

Since the two dimensional field theory is conformally invariant, if 0(t, A^) is a solution, 
then 0(t, 1/A^) is also a solution. This is clearly true for the solution in Eq. (26). Since 
A^ is a coordinate in our description, this connects up the infrared behaviour of the field 
theory with the ultraviolet behaviour. 

More directly, we note that the symmetries of the theory enhance significantly near 
the A^ = hypersurface. Conformal invariance, similar to the one found above, occurs 
in the gravitational sector as well. Defining q = — ^ by dq = —dN/N{Na), we see that 
A^ ~ exp(— Kg) near the horizon, where Na ~ k. The space part of the metric in Eq.(15) 
becomes, near the horizon dl"^ = N'^{dq^ + e'^'^'^dL\) which is conformal to the metric of 
the anti-De Sitter (AdS) space. The horizon becomes the g — )■ oo surface of the AdS 
space. These results hold in any dimension. There is a strong indication that most of the 
results related to horizons will arise from the enhanced symmetry of the theory near the 
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= surface (see e.g. [53, 54, 55, 56] and references cited therein). One can construct 
the metric in the bulk by a Taylor series expansion, from the form of the metric near 
the horizon, along the lines of Exercise 1 (page 290) of [57] to demonstrate the enhanced 
symmetry. These results arise because, algebraically, — ?■ makes certain terms in the 
diffeomorphisms vanish and increases the symmetry. This fact will prove to be useful 
in Sec. 5.2. 



2. 6. Three specific examples of horizons 

For the sake of reference, we briefly describe three specific solutions to Einstein's 
equations with horizons having a metric of the form in Eq. (8), viz. the Rindler, 
Schwarzschild and de Sitter spacetimes. In each of these cases, the metric can be 
expressed in the form of Eq. (8) with different forms of /(/) given in the Table 1. All 
these cases have only one horizon at some surface / = Ih and the surface gravity n is 
well defined. (We have relaxed the condition that the horizon occurs at / = 0; hence k 
is defined as (1/2)|/'| evaluated at the location of the horizon, / = l^-) The coordinate 
transformations relevant for analytic extension of these three spacetimes are also given 
in Table 1. The coordinates {T,X) are well behaved near the horizon while the original 
coordinate system (t, I) is singular at the horizon. Figure 1 describes all the three cases 
of horizons which we are interested in, with suitable definition for the coordinates. 

The horizons with the above features arise in Einstein's theory as well as in more 
general theories of gravity. While the detailed properties of the spacetimes in which 
these horizons occur are widely different, there are some key features shared by all the 
horizons we are interested in, which is worth summarizing: 

In all these cases, there exists a Killing vector field which is timelike in part 
of the manifold with the components = (1,0,0,0) in the Schwarzschild-type static 
coordinates. The norm of this field ^"'C,a vanishes on the horizon that acts as a bifurcation 
surface "H. Hence, the points of are fixed points of the Killing field. Further the surface 



Metric 



Rindler 



Schwarzschild 



De Sitter 



m 

^=l\f{lH)\ 

kX 
kT 



2kI 

j^ln kI 
'2kI cosh nt 
' IkI sinh Kt 



[1 



2n\ 



4M 

/ + 2Mln [^-1] 



r I 



1] 
1] 



1/2 



1/2 



cosh [^] 



(1 - HH , 

2H \1+HU 

imf'-oshHt 

(l^)^/%inhm 



Table 1. Properties of Rindler, Schwarzschild and De Sitter metrics 
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gravity k of the horizon can be defined using the 'acceleration' of the Kilhng vector by: 

^^VfeC = (28) 
When defined in this manner, the value of n depends on the normalization chosen for 
(If we rescale C,"" — )■ fJ'C,"', the surface gravity also scales as k — )■ fin). Very often, 
however, we will be interested in the combination C,°'/k,, which is invariant under this 
scaling. 

In these spacetimes, there exists a spacelike hypersurface S which includes "H and 
is divided by "H into two pieces and E^, the intersection of which is in fact T-L. In 
the case of black hole manifold, for example, E is the T = surface, E/j and E^, are 
parts of it in the right and left wedges and "H corresponds to the / = 2M surface. The 
topology of Eij and "H depends on the details of the spacetime but "H is assumed to have 
a non-zero surface gravity. Given this structure it is possible to generalize most of the 
results we will be discussing in the later sections. 

Finally, to conclude this section, we shall summarize a series of geometrical facts 
related to the Rindler frame and the Rindler horizons which will turn out to be useful 
in our future discussions. Though we will present the results in the context of a 2- 
dimensional Rindler spacetime, most of the ideas have a very natural generalization to 
other bifurcation horizons. We begin with the metric for the Rindler spacetime expressed 
in different sets of coordinates: 

dp 

ds^ = -dT^ + dX^ = -dUdV = -e^^^'-^'^dudv = -2kI dt^ + — (29) 

The coordinate transformations relating these have been discussed earlier. In particular, 
note that the null coordinates in the two frames are related hy U = T — X = 
— /t~^e~''", V = T + X = K~^e'^'" . We will now introduce several closely related vectors 
and their properties. 

(i) Let k"" be a future directed null vector with components proportional to (1,1) in 
the inertial frame. The corresponding affinely parameterized null curve can be taken to 
be x°' = kX{1, 1)j with X being the affine parameter and subscript / (or R) indicates 
the components in the inertial (or Rindler) frame. 

(ii) We also have the natural Killing vector corresponding to translations in the 
Rindler time coordinate. This vector has the components, = (l,0)/j = k,{X,T)j and 
^""^a = —'2,kI = —N"^. This shows that the bifurcation horizon T-i is at the location where 
i"'ia = 0. The "acceleration" of this Killing vector is given by a* = .^''V;,^* = k^(T, X)/ 
and hence, on the horizon, a* = nC,^ consistent with Eq. (28). It is also easy to see that, 
on the horizon, — nXk^ with k"" = (1, 1)/. 

(iii) Another natural vector which arises in the Rindler frame is the four-velocity 
u"" of observers, moving along the orbits of the Killing vector C,"'. On the horizon "H, this 
four-velocity has the limiting behaviour Nu"" — nXk"". 

(iv) Lastly, we introduce the unit normal to / = constant surface, which also has 
the limiting behaviour Nra — )■ nXk"' when we approach the horizon. It therefore follows 
that Nu^,Nr^,a^ and all tend to vectors proportional to fc* on the horizon. These 
facts will prove to be useful in our later discussions. 
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3. Thermodynamics of horizon: A first look 

We shall now provide a general argument which associates a non-zero temperature with 
a bifurcation horizon. This argument, originally due to T.D. Lee, ([58]; also see ref. [59]) 
is quite powerful and elegant and applies to all the horizons which we will be interested 
in. It uses techniques from path integral approach to quantum field theory, which we 
shall first review briefly. 

3.1. Review of Path Integral approach 

It is known in quantum mechanics that the net probability amplitude K{2] 1) for the 
particle to go from the event Vi to the event V2 is obtained by adding up the amplitudes 
for all the paths connecting the events: 

K{V2; Vi) = K{t2, ga; h, gi) = J2 exp[^^(path)] (30) 

paths 

where ^(path) is the action evaluated for a given path connecting the end points Vi 
and 7^2 • The addition of the amplitudes allows for the quantum mechanical interference 
between the paths. The quantity K{t2, q2', ti, qi) contains the full dynamical information 
about the quantum mechanical system. Given K{t2, q2] ti, qi) and the initial amplitude 
ijj{ti, qi) for the particle to be found at gi, we can compute the wave function ip(t, q) at 
any later time by the usual rules for combining the amplitudes: 

V'(i, 9) = y dqiK{t, q; ti, qi)^{ti, qi) (31) 

The above expressions continue to hold even when we deal with several degrees of 
freedom qi,q2,... which may still be collectively denoted as g = {qi}; it is understood 
that the integral in Eq. (31) has to be performed over all the degrees of freedom. 

To obtain the corresponding results in field theory, one needs to go from a discrete 
set of degrees of freedom (labeled by i = 1,2,....) to a continuum of variables denoting 
the coordinates x in a spacelike hypersurface. In this case the dynamical variable at 
time t = ti is the field configuration g(x). (For every value of x we have one degree of 
freedom.) The integral in Eq. (31) now becomes a functional integral over the initial 
field configuration and Eq. (31) becomes 

ij{t, g(x)) = j VqiKit, g(x); ti, gi(x))V^(ti, gi(x)) (32) 

We shall, however, not bother to indicate this difference between field theory and point 
quantum mechanics and will continue to work with latter since the generalizations will 
be quite obvious by context. 

We will next obtain a relation between the ground state wave function of the system 
and the path integral kernel which will prove to be useful. In the conventional approach 
to quantum mechanics, using the Heisenberg picture, we will describe the system in 
terms of the position and momentum operators g and p. Let \q,t) be the eigenstate 
of the operator q{t) with eigenvalue g. The kernel — which represents the probability 
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amplitude for a particle to propagate from to (t2,q2) — can be expressed, in a 

more conventional notation, as the matrix element: 

K(t2,q2;ti,qi) = {q2,t2\ti,qi) = {q2,0\exp[-iH(t2 -ti)]\0,qi). (33) 

where H is the time-independent Hamiltonian describing the system. This relation 
allows one to represent the kernel in terms of the energy eigenstates of the system. We 
have 

K{T, ga; 0, gi) = (ga, 0| exp -iHT\0, qi) 

= '^{q2\En){En\exp-iHT\Em){Em\qi) 

n,m 

= J2MQ2)rn{qi)eM-^Er,T) (34) 

n 

where il^n{q) = {q\En) is the n-th energy eigenfunction of the system under consideration. 
Equation (34) allows one to express the kernel in terms of the eigenfunctions of the 
Hamiltonian. For any Hamiltonian which is bounded from below it is convenient to add 
a constant to the Hamiltonian so that the ground state — corresponding to the n = 
term in the above expression — has zero energy. We shall assume that this is done. 
Next, we will analytically continue the expression in Eq. (34) to imaginary values of T 
by writing iT = Te- The Euclidean kernel obtained from Eq. (34) has the form 

Ke{Te, q2] 0, gi) = ^ V'n('?2)V'n('?i) exp(-KTB) (35) 

n 

Suppose we now set gi = g, ga = in the above expression and take the limit Tg; — )■ oo. 
In the large time limit, the exponential will suppress all the terms in the sum except 
the one with En = ^ which is the ground state for which the wave function is real. We, 
therefore, obtain the result 

lim K{T, 0; 0, g) ^ ^o(0)^o(gi) « ^o(g) (36) 

T->-oo 

Hence the ground state wave function can be obtained by analytically continuing the 
kernel into imaginary time and taking a suitable limit. The proportionality constant in 
Eq. (36) is irrelevant since it can always be obtained by normalizing the wave function 
ipo{q). Hence we have 

i){q) oc K (oo, 0; 0, g) = K (0, g; -cx), 0) (37) 

where, in the arguments of the first one refers to Euclidean time and the second one 
refers to the dynamical variable. The last equality is obtained by noting that in Eq. (36) 
we can take the limit T oo either by (ta — > oo, ti = 0) or by (ta = 0, ti — — oo). This 
result holds for any closed system with bounded Hamiltonian. Expressing the kernel as 
a path integral we can write this result in the form 

JTE=0,q 

This formula is also valid in field theory if g is replaced by the field configuration g(x). 
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The analytic continuation to imaginary values of time also has close mathematical 
connections with the description of systems in thermal bath. To see this, consider the 
mean value of some observable 0{q) of a quantum mechanical system. If the system 
is in an energy eigenstate described by the wave function ipniq), then the expectation 
value of 0{q) can be obtained by integrating 0{q)\ipn{q)\'^ over q. If the system is in 
a thermal bath at temperature (3~^, described by a canonical ensemble, then the mean 
value has to be computed by averaging over all the energy eigenstates as well with a 
weightage exp(— In this case, the mean value can be expressed as 

= / dqMq)0{q)rn{q)e-^^" f dqp{q,q)0{q) (39) 

where Z is the partition function and we have defined a density matrix p{q, q') by 

p{q,q')^Y.'^-{q)rniq')e-^''" (40) 

n 

in terms of which we can rewrite Eq. (39) as 

where the trace operation involves setting q = q' and integrating over q. This standard 
result shows how p(g, q') contains information about both thermal and quantum 
mechanical averaging. Comparing Eq. (40) with Eq. (34) we find that the density 
matrix can be immediately obtained from the Euclidean kernel by: 

p(g,g') = ^i^(/3,g;0,g') (42) 
with the Euclidean time acting as inverse temperature. 



3.2. Horizon temperature from a path integral 

We shall now consider the quantum field theory in a spacetime with a horizon which 
can be described in two different coordinate systems. The first one (T, X, x_l) is a 
global coordinate system which covers the entire spacetime manifold (which could be 
the inertial Cartesian coordinate system in flat spacetime or the Kruskal-like coordinate 
system in the case of spherically symmetric metrics with horizon). The second one 
(t, x,xx) covers the four different quadrants of the spacetime and is related to the first 
set by a set of transformations similar to Eq. (9). We shall now show that the global 
vacuum state defined on the surface T = appears as a thermal state to observers 
confined to the right wedge TZ with a temperature = (27r/fi;). 

On analytic continuation to imaginary time, the two sets of coordinates behave as 
shown in Fig. 2. The key new feature is that tE becomes an angular coordinate having a 
periodicity {2tx/ hi). While the evolution in (effected by the inertial Hamiltonian Hi) 
will take the field configuration from Tj; = to T^; — )■ oo, the same time evolution gets 
mapped in terms oHe into evolving the "angular" coordinate tE from to 27r/K and is 
effected by the Rindler Hamiltonian H^i. (This should be clear from Fig. 2.) It is obvious 
that the entire upper half-plane T > is covered in two completely different ways in 




terms of the evolution in Te compared to evolution in Ie- In (Te,X) coordinates, we 
vary X in the range (—00, 00) for each Te and vary Te in the range (0, 00). In {tE,x) 
coordinates, we vary x in the range (0, 00) for each tE and vary Ie in the range (0, vt/k). 
This fact allows us to prove that 

(vac|0z., <Pr) cx (0L|e--^«/^|0B) (43) 

as we shall see below. 

To provide a simple proof of Eq. (43), let us consider the ground state wave 
functional (vac|0i,0fl) in the extended spacetime expressed as a path integral. From 
Eq. (38) we know that the ground state wave functional can be represented as a 
Euclidean path integral of the form 

(vac|g) oc / V^e-^ (44) 

JTE=0;^=q 

where Te = iT is the Euclidean time coordinate and we have denoted the field 
configuration on the T = hypersurface by (^(x). But we know that this field 
configuration can also be specified uniquely by specifying 0l(x) = g(x) with X < 
and 0_r(x) = g(x) with X > 0. Hence we can write the above result in terms of 0^ and 

(t)R as 

;.Tb=cx);0={O,O) 

{ys.c\<PlAr) oc / P0e-^ (45) 

JTE=0;<t>=(<t>L,<t>R) 

From Fig. 2 it is obvious that this path integral could also be evaluated in the polar 
coordinates by varying the angle 9 = kIe from to tt. When ^ = the field configuration 
corresponds to = 0^ and when 9 = tt the field configuration corresponds to = 0^. 
Therefore Eq. (45) can also be expressed as 

(vac|0L,0R) oc / P0e-^ (46) 

J KtE=0;4'=4>R 
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But in the Heisenberg picture, 'rotating' from nt^ = to nt^ = ir is a. time evolution 
governed by the Rindler Hamiltonian Hr. So the path integral Eq. (46) can be 
represented as a matrix element of the Rindler Hamiltonian Hr giving us the result: 

/•KtE=iT;<l>=(j>L 

(vac|0z.,0^) oc / VcPe"^ = (0L|e-(-/'')^«|0^) (47) 

J K,tE=0;<j>=4>[i 

proving Eq. (43). 

If we denote the proportionality constant in Eq. (43) by C, then the normalization 
condition 

1 = y" 'D(f)L'D(j)R\{vac\(j)L,(l)R)f = j I'0L'P0i?. (vac|0L,0ij)(0L,0R|vac) 

(48) 

fixes the proportionality constant C, allowing us to write Eq. (43) in the form: 

<""^''^^''^^^= [Tr(e-^-^A)]^/^ ^ ^ 

From this result, we can compute the density matrix for observations confined to the 
Rindler wedge 71 by tracing out the field configuration 0l on the left wedge. We get: 



p{4>R, (P'r) = j I'0L(vac|0i, 0'j?|vac) 



7-^(g-27rHfl/K-) ^^^^ 

Thus, tracing over the field configuration 0l in the region behind the horizon leads to 
a thermal density matrix p oc exp[— (27r/K)if/j] for the observables in TZ. 

The main ingredients which have gone into this result are the following, (i) The 
singular behaviour of the (t, x) coordinate system near x = divides the T = 
hypersurface into two separate regions, (ii) In terms of real (t, x) coordinates, it is 
not possible to distinguish between the points (T, X) and (— T, — X) but the complex 
transformation t ^ t ± in maps the point {T,X) to the point (— T, —X). That is, a 
rotation in the complex plane (Re t, Im t) encodes the information contained in the full 
T = plane. 

In fact, one can obtain the expression for the density matrix directly from path 
integrals along the following lines. We begin with the standard relation in Eq. (37) 
which gives 

(vac|g) = K{oo, 0; 0, q) = K{0, q; -cx), 0) = (vac|0L, 0r) (51) 

where, in the arguments of K, the first one refers to Euclidean time and the second one 
refers to the dynamical variable. The density matrix used by the observer in the right 
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(52) 



where we explicitly decomposed g(x) into the set ^/.(x), ^^^(x) everywhere. The 
expression in the right hand side evolves the system from Te = —oo to Te = +00 
with some specific restrictions on the field configuration on the Te = hypersurface. 
Since we are integrating out all the field configurations (pi, it follows that there is no 
restriction on the field along X < 0. To handle the field configurations on X > 0, 
we can proceed as follows. We first note that = is the same as t^; = when 
X > 0. Instead of considering = t^; = 0, let us consider an infinitesimally displaced 
hypersurface t^; = e in one of the kernels and tE = (27r//t) — e in the second kernel. That 
is, instead of specifying 0_r(x) and 0^(x) at = we will specify 0_r(x) at tE = e and 
4>'p>{^) at tE = {2tt / k) — e. It is then obvious that the result of the integral in Eq. (52) 
is the propagation kernel that propagates 0_r(x) at t^; = e to 0^(x) at tE = {27i/k) — e 
in the Rindler time. Taking the e — )■ limit is now trivial and we get 



where Hji is the Rindler Hamiltonian. Comparing the first and last expressions we find 
that the operator corresponding to the density matrix is just p = exp[—{27i/K)Hii]. 
(This is an unnormalized density matrix since we have not bothered to normalize the 
wavef unctions.) 

3.3. Complex time and the region beyond the horizon 

There are two points that need to be stressed regarding the above derivation. First, the 
light cones described by X"^ — = get mapped to the origin of the Euclidean sector 
through the equation X"^ + T| = 0. Consequently the quadrants J-" and V disappear 
from the Euclidean sector — or rather, they collapse into the origin. This implies that 
the region beyond the horizon is not covered by the Euclidean coordinates. The second 
point is that even though the horizon collapses to a single point in the Euclidean sector, 
the Euclidean Rindler time tE contains information about the left quadrant C To see 
this, we only have to compare Eq. (12) taken with the positive sign and Eq. (17). When 
t ranges from —00 to +00, the coordinate X = xcosh/tt remains positive. On the other 
hand, when the corresponding Euclidean time tE varies from to 27r/K, the coordinate 
X = X cos KtE covers both the right wedge TZ and the left wedge C. Therefore the range 
of Euclidean Rindler time from tE = t^/2k to tE = —ti/2k covers the region beyond the 
horizon [60]. It is because of this peculiar feature (which, of course, is closely related 
to periodicity in the imaginary time) that we can obtain the thermal effect due to 




(53) 



i^((27r/«:),0^;O,0R) = ((/)'«| exp[-(27r/«:)i/R] 
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horizon from the Euchdean approach. The conclusions are strengthened by a few other 
considerations which are worth mentioning briefly. 

To begin with, note that similar results arise in a more general context for 
any system described by a wave function \Ef(t,/;£') = exp[iA(t,l; E)] in the WKB 
approximation [61]. The dependence of the quantum mechanical probability P{E) = 
|\E'P on the energy E can be quantified in terms of the derivative 

d\nP d , (dA\ , , 

-— 2(ImA) = -2Im — (54) 



dE dE ^ ' \dE ^ 

in which the dependence on (t, /) is suppressed. Under normal circumstances, action will 
be real in the leading order approximation and the imaginary part will vanish. (One 
well known counter-example is in the case of tunneling in which the action acquires 
an imaginary part; Eq. (54) then correctly describes the dependence of tunneling 
probability on the energy.) For any Hamiltonian system, the quantity [dAjdE] can 
be set to a constant to to determine the trajectory of the system: [dA/dE) = — to- 
Once the trajectory is known, this equation determines to a function of E [as well as 
(t, /)]. Hence we can write 
dlnP 

^^2lm[tom (55) 

From the trajectory in Eq. (18) which is valid near the horizon, we note that to{E) can 
pick up an imaginary part if the trajectory of the system crosses the horizon. In fact, 
since nt ^ nt — in changes X to —X [see Eqs. (5,6,9)], the imaginary part is given by 
(— tt/k) leading to [din P/dE) = —27i/k. Integrating, we find that the probability for 
the trajectory of any system to cross the horizon, with the energy E, will be given by 
the Boltzmann factor 

2n 



P{E) (X exp 



-E 

K 



Poexp[-/3E] (56) 



with temperature T = nj^ix. (For special cases of this general result see [62] and 
references cited therein.) 

It is also interesting to examine how these results relate to the more formal approach 
to quantum field theory. The relation between quantum field theories in two sets of 
coordinates (t, x) and (T, X), related by Eq. (9), with the metric being static in the (t, x) 
coordinates can be described as follows: Static nature suggests a natural decomposition 
of wave modes as 

0(t, x) = I rfu;[a^/.(x)e-'^* + ai/:(x)e^-*] (57) 

in (t,x) coordinates. These modes, however, behave badly (as x^'"^l'^\ see Eq. (26)) 
near the horizon since the metric is singular near the horizon in these coordinates. We 
could, however, expand 0(t, x) in terms of some other set of modes F^(t,x) which are 
well behaved at the horizon. This could, for example, be done by solving the wave 
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equation in (T, X) coordinates and rewriting the solution in terms of (t,x). This gives 
an alternative expansion for the field: 

x) = y diy[A,F,{t, x) + AlF:{t, x)] (58) 

Both these sets of creation and annihilation operators define two different vacuum states 
O'ui\0)a = 0,^4^10)^ = 0. The modes F^(t,x) will contain both positive and negative 
frequency components with respect to t while the modes /^(x)e~*"* are pure positive 
frequency components. The positive and negative frequency components of Fi,(t,x) can 
be extracted through the Fourier transforms: 

/CO /"OO 
dt e''^'F,{t,Xf); (3^,= / dt e-''''F,{t,^j) (59) 
oo J — oo 

where xj is some convenient fiducial location far away from the horizon. One can think 
of looji/P and as similar to unnormalized transmission and reflection coefficients. 

(They are very closely related to the Bogoliubov coefficients usually used to relate 
two sets of creation and annihilation operators.) The a— particles in the \0)a state 
is determined by the quantity l/S^jiy/ctwiyp- If the particles are uncorrelated, then the 
standard relation iV|Q;p = (A^ + l)|/9p between absorption and emission leads to the 
flux of out-going particles: 

N = (60) 
1-1/3 /a P ^ ' 

If the F modes are chosen to be regular near the horizon, varying as exp{—iQU) etc., 
then Eq. (9) shows that F^{t,:}Cf) oc exp{—iQqe~'^^) etc. The integrals in Eq. (59) 
again reduces to the Fourier transform of an exponentially redshifted wave and we get 
l/^uju/dwul'^ = e"^'^ and Eq. (60) leads to the Planck spectrum. This is the quantum 
mechanical version of Eq. (20) and Eq. (23). 

Finally, one can relate the above result to the analyticity properties of wave modes 
of a scalar field in the U = (T — X) coordinates and u = (t — x) coordinates. Since the 
positive frequency mode solution to the wave equation in the (T, X) coordinates has the 
form exp{—iQU) (with 1] > 0) and is analytic in the lower half of complex U plane, any 
arbitrary superposition of such modes with different (positive) values of Q will also be 
analytic in the lower half of complex U plane. Conversely, if we construct a mode which 
is analytic in the lower half of complex U plane, it can be expressed as a superposition of 
purely positive frequency modes [29]. From the transformations in Eq. (9), we find that 
the positive frequency wave mode near the horizon, = exp{—iuu) can be expressed 
as (p (X U'^^/'^ for U < If we interpret this mode as oc (?7 — ief'^l'^ then, this mode 
is analytic throughout the lower half of complex V plane. We can then interpret the 
mode as 

f e[^("/")'"^] (for f/ > 0) 

|eW«e[(W«)in|f/|) (forf/<0) 

This interpretation of ln(— t/) as In |f/| — ivr = km — ivr = /tt — ^ — iyr is consistent with 
the procedure of replacing nt ^ nt — iir to go from X > to X < 0. This is precisely 
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what happens in the Euchdean continuation. The factor e'^"/'' in the second hne of the 
Eq. (61) leads to the thermal effects in the conventional picture. 

3.4- Hawking radiation from black holes 

The description in the previous sections shows that the vacuum state defined in a 
coordinate system — which covers the full manifold — appears as a thermal state 
to an observer who is confined to part of the manifold partitioned by a horizon. This 
result (and the analysis) will hold for any static spacetime with a bifurcation horizon, 
like the Schwarzschild spacetime, de Sitter spacetime etc. All these cases describe a 
situation in thermal equilibrium at a temperature T = k/2tt (where k is the surface 
gravity of the horizon) as far as an observer confined to the region R is concerned. 

A completely different phenomenon arises in the case of a dynamical situation like, 
for example, the collapse of a spherically symmetric massive body to form a black hole. 
In this case, time reversal invariance is explicitly broken. The study of a quantum field 
theory in such a context shows that, at late times, there will be a flux of radiation 
flowing towards the future null infinity with a Planckian spectrum corresponding to a 
temperature K/2n. This process is called black hole evaporation. 

This result is conceptually different from associating a temperature with the 
horizon. In the case of a Rindler spacetime, for example, there is no steady flux of 
radiation propagating towards future null infinity even though an observer confined 
to the region TZ will interpret the expectation values of operators as thermal averages 
corresponding to a temperature k/27i. This corresponds to a situation which is time 
reversal invariant characterized by thermal equilibrium. The black hole evaporation, in 
contrast, is an irreversible process. 

We shall now work out the corresponding result for a scalar field in the time 
dependent metric generated by collapsing matter. The scalar field can be decomposed 
into positive and negative frequency modes in the usual manner. We choose these modes 
in such a way that at early times they correspond to a vacuum state. In the presence of 
collapsing matter, these modes evolve at late times to those with exponential redshift, 
thereby leading to thermal behaviour. 

To do this, we need an explicit model for the collapsing matter. Since only the 
exponential redshift of the modes at late times is relevant as far as the thermal spectrum 
is concerned, the result should be independent of the detailed nature of the collapsing 
matter. So we shall choose a simple model for the formation of the black hole, based 
on a spherical shell of mass M that collapses under its own weight. The metric inside 
the shell will be fiat while the one outside will be Schwarzschild. Further, the angular 
coordinates do not play a significant role in this analysis, allowing us to work in the 
2-dimensional (t, r) subspace. 

The line element outside and inside the collapsing, spherically symmetric, collapsing 
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shell is taken to be 
ds^ = 



-C{r) dudv 
-dUdV 



2M 



dudv 



(exterior) 
(interior) 



where 



u 



t-^ + m; v = t + ^-m; e 



dr C-^- 



U = T-r + Ro, 



V = T + r-Ro 



(62) 



(63) 
(64) 



The -Ro and Rq are constants related in the same manner as r and C,- 

Let us assume that, for r < 0, matter was at rest with its surface aX r = Rq and 
for r > 0, it collapses inward along the trajectory r = R{t). The coordinates have been 
chosen so that at the onset of collapse (r = t = 0) we have u = U = v = V = Q aX the 
surface. Let the coordinate transformations between the interior and exterior be given 
by the functional forms U = f{u) and v = h{V). Matching the geometry along the 
trajectory r = R{t), requires the condition 

^ du dV / „ / X X /X 

Using Eq. (63) and Eq. (64) along the trajectory, this equation can be simplified to give 

2 7S2 



dt\- 
d^J 



— 1(1 



R' 



(66) 



where R denotes dR/dr and f/, V and C are evaluated along r = R{t). Using this again 
in the definition of u, v etc., it is easy to show that 



dU_ 

du 

dv 
dV 



dl 
du 

dh 
dV 



l-R)Ci C{l-R 



C(l - R" 



R' 



1/2 



C{l + R) 



R' 



1/2 



R 



(67) 



(68) 



Since i? < for the collapsing shell, we should take (i?^)^/^ = — i?. 

We now introduce a massless scalar field in this spacetime which satisfies the 
equation 00 = 0. As the modes of the scalar field propagate inwards they will reach 
r = and re-emerge as out-going modes. In the (t, r) plane, this requires reflection of 
the modes on the r = line, which corresponds io V = U — 2Rq. Since the modes 
vanish for r < 0, continuity requires = at r = 0. The solutions to the 2-dimensional 
wave equations = which (i) vanish on the line V = U — 2Rq and (ii) reduce to 
standard exponential form in the remote past, can be determined by noting that, along 
r = we have 



V = h{V) = h[U - 2Ro] = h[f{u) - 2Ro] 



(69) 
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where the square bracket denotes functional dependence of h on its argument. Hence 
the solution is 

Given the trajectory -R(r), one can integrate Eq. (67) to obtain j{u) and use Eq. (70) 
to completely solve the problem. This will describe time-dependent particle production 
from some collapsing matter distribution and — in general — the results will depend 
on the details of the collapse [35, 63]. 

The analysis, however, simplifies considerably and a universal character emerges if 
the collapse proceeds to form a horizon on which C — ?■ 0. Near C = 0, equations (67) 
and (68) simplify to 

rf^^^^^^^' '^^~m~ ^^^^ 

where we have used the fact that (i?^)^/^ = —R for the collapsing solution. Further, 
near C = 0, we can expand -R(r) as -R(t) = v{jh — t) + 0\{rh — t)^] where R = Rh 
at the horizon and v = —R{Th). We have denoted by Th the time at which the shell 
crosses the horizon. Integrating Eq. (71) we get 

Ku ^ — In \ U + Rh — Rq — Th\ + const (72) 

where k = {l/2){dC / dr)^^^ is the surface gravity and 

V ^ constant -V{1 + v)/2u (73) 

It is now clear that: (i) The relation between v and V is linear and hence holds no 
surprises, (ii) The relation between U and u, which can be written as f/ oc exp(— km) 
signifies the exponential redshift we have alluded to several times. The late time 
behaviour of out-going modes can now be determined using Eq. (72) and Eq. (73) 
in Eq. (70). We get: 

$ ^ (e-*^'' - exp [iuj [ce"'^" + d])) (74) 

where c, d are constants. This mode with exponential redshift, can be expressed in terms 
of the modes exp(=l=i/M) as 

/■°° dy 

Determining a^ji/, Puv by Fourier transforming this relation, we get 

a^, = = ^— e-^I^^Y{-ivlK)- = e--'l^al, (76) 

'l.'KK^VbJ \ UJC J 

Note that these quantities do depend on c, d etc; but the modulus 

'^"'^'^ " 2;;:[exp(27rz//K)-l)] ^^^"^ 

is independent of these factors. (The mathematics is essentially the same as in Eqs.(22), 
(23)). This shows that the vacuum state at early times will be interpreted as containing 
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a thermal spectrum of particles at late times with temperature T = h/2ti. In the case 
of a black hole, k = (1/4M) and the temperature turns out to be T = (l/SvrM). 

Our result implies that when spherically symmetric configuration of matter 
collapses to form a black hole, observers at large distances will receive a thermal radiation 
of particles from the black hole at late times. (It is possible to prove this more formally 
by considering the expectation values of the energy-momentum tensor of the scalar field; 
this will demonstrate the flux of energy to large distances). It seems natural to assume 
that the source of this energy radiated to infinity is the mass of the collapsing structure. 
Strictly speaking, this is an extrapolation from our result because it involves changes 
in the background metric — which was parametrized by M - due to the effect of the 
radiation while our original result was based on a test scalar field in a fixed background 
metric. We shall, nevertheless, make this assumption and explore its consequences. 
Given the temperature of the black hole T{E) = I/SttM as a function of the energy 
E = M, we can integrate the expression dS = dE/T{E) to define an 'entropy' S{E) for 
the black hole: 



where Ahor = 47r(2M)^ is the area of the r = 2M, t = constant surface and Lp = Gh/c^ 
is the square of the Planck length. This result shows that the entropy obtained by 
integrating the T{E) is proportional to the area of the horizon. 

This result connects up with several classical features of black holes. We know 
that the area of the horizon does not decrease during classical processes involving 
the black holes which suggests an analogy between horizon area and entropy. The 
Eq. (2) for classical processes involving black holes now acquires a direct thermodynamic 
interpretation. The factor {k/2'k) and {A/A) can indeed be identified with physical 
temperature and entropy. Note that classical analysis can only identify these quantities 
up to a multiplicative factor. On the other hand, the analysis of quantum fields in the 
Schwarzschild metric allows us to determine the temperature and entropy uniquely and 
the entropy turns out to be one-quarter of the area of the horizon expressed in Planck 
units. § 

3. 5. Horizon entropy in generalized theories of gravity 

There has been considerable amount of work in analyzing the nature of horizon 
thermodynamics in theories different from general relativity. In a wide class of such 
theories, one does get solutions with horizons and one can associate a temperature and 
entropy with them. While the temperature can be identified from the periodicity of 

§ There is an ambiguity in the overall additive eonstant to entropy whieh is settled in Eq. (78) by 
assuming that 5* = for M = 0. This might appear reasonable but recall that T — > oo when M — >■ 0; 
flat spacetime, treated as the Af — > limit of Schwarzschild metric, has infinite temperature rather 
than zero temperature. Hence, it is worth emphasizing that choosing S ~ for M = is a specific 




(78) 



assumption. 
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Euclidean time, determining the correct form of entropy is more non-trivial. We shall 
now briefly describe how these results arise in a class of theories which are natural 
generalizations of Einstein gravity. 

Consider a theory for gravity described by the metric gab coupled to matter. We 
will take the action describing such a theory in D— dimensions to be 

A = J d^x,/^ [L{Rf,, g-') + Lrnauig"', Qa)] (79) 

where L is any scalar built from metric and curvature and Lmatt is the matter Lagrangian 
depending on the metric and some matter variables g^. (We have assumed that L does 
not involve derivatives of curvature tensor, to simplify the discussion; a more general 
structure is explored e.g., in ref. [64]) Varying g"-^ in Eq. (79) we get 5(Lmattv^— ^) = 
-{l/2)^Tab5g''^ and 

8{L^) = {GabSg'''' + V.^^;'^) . (80) 

The variation of the gravitational Lagrangian density generically leads to a surface term 
which is expressed by the Va(5f") term. Ignoring this term for the moment (we will 
comment on this later) we get equations of motion (see e.g. Refs. [41, 65]) to be 
2Qab = Tab where the exphcit form of Qab is 

Gab = Pa^'-Rbcde ' \Lgab " IV^V'Paaib ^ T^ab - 2V^V'Pacdb (81) 

where 

pabcd ^ (82) 
dRabcd 

(Our notation is based on the fact that Qab = Gab, T^ab = Rab in Einstein's gravity.) For 
any Lagrangian L, the functional derivative Qab satisfies the generalized off-shell Bianchi 
identity: VaQ"^ = 0. 

Many such models have been investigated in the literature and most of these models 
have black hole solutions. Whenever the black hole metric can be approximated by a 
Rindler metric near the horizon, it is possible to associate a temperature with the 
horizon, using the procedures described earlier, e.g. in Sec. (3.2). Associating the 
entropy is more nontrivial and we shall now indicate how this is usually done. (For 
a rigorous proof, which we shall not provide, see e.g., ref [11]. We will, however, provide 
an alternative route to this result in Section 5.8.) 

In any generally covariant theory, the infinitesimal coordinate transformations 
x°- x°- + C,"" lead to conservation of a Noether current which depends on C,"". To derive 
the expression for the Noether current, let us consider the variations in 6 gab which arise 
through the diffeomorphism x"^ — ?■ + In this case, 6{Ly/—g) = —^/—g'Va{L^°'), 
with Sg'''' = (V"^^ + V^^'^). Substituting these in Eq. (80) and using Va^"* = 0, we 
obtain the conservation law VaJ"" = 0, for the current, 

ja ^ (^2Q''''^b + LC + ^gt;") = 27^'^''^6 + 5^1;" (83) 
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where 6^v°- represents the boundary term which arises for the specific variation of the 
metric in the form Sg""^ = (V".^* + V^.^"). Quite generally, the boundary term can be 
expressed as [65], 

5v^ = ^-a^^"^^5g,^ + l(3f'^5Ti (84) 

where, we have used the notation Q*-' = Q^-' + Q-'*. The coefficient p"-^'^'^ arises from the 
derivative of Lgrav with respect to R"-^'^'^ and hence possess all the algebraic symmetries 
of curvature tensor. In the special case of diffeomorphisms, x"" — )■ + the variation 
" is given by Eq. (84) with: 

6gab = -^ia^b); 5r^e=-^V(,V.)e' + ^/2?,e)™r (85) 

Using the above expressions in Eq. (83), it is possible to write an explicit expression 
for the current J" for any diffeomorphism invariant theory. It is also convenient to 
introduce an anti-symmetric tensor J"'' by J" = V^J'^*. For the general class of theories 
we are considering, the J"'^ and J" can be expressed in the form 

jab ^ ^pabcdy^^^ _ {y^pabcd^^ (gg) 

r = -2Vb {P'"^''^ + P'^'*'^) Vc^d + 2P'^''"'^VfeVc^d - A^dVbVcP''''"^ (87) 

where Pabcd = {dL/dW'^^'^). (The expression for J^, J ah are not unique. This ambiguity 
has been extensively discussed in the literature but for our purpose we will use the J*^ 
defined as above.) 

We shall see that, for most of our discussion, we will not require the explicit form 
of d^v"" except for one easily proved result: d^v"" = when C,"" is a Killing vector and 
satisfies the conditions 

V(a6) = 0; VaVbCc = Rcbad^'' (88) 

The expression for Noether current simplifies considerably when satisfies Eq. (88) 
and is given by 

ja ^ + LC) = 2Tr^ih (89) 

The integral of J'^ over a spacelike surface defines the conserved Noether charge, M . 

To obtain a relation between the horizon entropy and Noether charge, we first note 
that on-shell i.e., when field equations hold {2Qab = Tab), we can write: 

ja ^ ^j^aj ^ gajj^^ ^, (gg^ 

Therefore, for any vector ka which satisfies ka^"" = 0, we get the result: 

(kar) = T'^'kaij. (91) 

The change in this quantity, when T""^ changes by a small amount dT""^ , will be 
5{kaJ'^) = kaC,jST"-^ . It is this relation which can be used to obtain an expression 
for horizon entropy in terms of the Noether charge. When some amount of matter with 
energy-momentum tensor crosses the horizon, the corresponding energy fiux can 
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be thought of as given by integral of kaC,jST"-^ over the horizon where is the Kilhng 
vector field corresponding to the bifurcation horizon and ka is a vector orthogonal to 
which can be taken as the normal to a timelike surface, infinitesimally away from the 
horizon. (Such a surface is sometimes called a 'stretched horizon' and is defined by the 
condition N = e where N is the lapse function with iV = representing the horizon.) 
In the [D — 1) dimensional integral over this surface, one coordinate is just time; since 
we are dealing with an approximately stationary situation, the time integral reduces 
to multiplication by the range of integration. Based on our discussion earlier we will 
assume that the time integration can be restricted to the range (0, /3) where /3 = 27r/fi; 
and K is the surface gravity of the horizon. (The justification for this requires a much 
more detailed mathematical analysis which we shall not get into.) Thus, on integrating 
5{kaJ°') over the horizon we get 

5 [ d^-^xVhikaJ") = [ d^~^xVhka^j6T^^ 

JH JH 

= f3 [ d^-^xVhka^jST"^ (92) 

JH 

where the integration over time has been replaced by a multiplication by /3 = (2-11 /n) 
assuming approximate stationarity of the expression. The integral over ^T"-' is the flux 
of energy 6E through the horizon so that PSE can be interpreted as the rate of change 
of the entropy associated with this energy flux. One can obtain, using these facts, an 
expression for entropy, given by 

jD-1v 70. _ P f ^D-2^r^ jab 



^Noether = ^ = j d^'-^T.^r = '- j d^'-'T.abr' (93) 

where d^~^T,a = d^^^xy/hka, the Noether charge is A/" and we have introduced the 
antisymmetric tensor J"^ by J" = V^J"^. In the final expression the integral is over any 
surface with {D — 2) dimension which is a spacelike cross-section of the Killing horizon 
on which the norm of vanishes. 

As an example, consider the special case of Einstein gravity in which Eq. (86) 
reduces to 

jab ^ _±_ /ya^6 _ ^b^a\ .g^N 

If be the timelike Killing vector in the spacetime describing a Schwarzschild black 
hole., we can compute the Noether charge Af as an integral of J""^ over any two surface 
which is a spacelike cross-section of the Killing horizon on which the norm of vanishes. 
The area element on the horizon can be taken to be dHab = {^a^b — h^a)\f^d^~'^x in 
Eq. (93) with la being an auxiliary vector field satisfying the condition Za^" = — 1. Then 
the integral in Eq. (93) reduces to 



S Noether — ~ TT" / \f^d^ '^^i^a^b)^''^"' 



g f V^d^-'x = \An (95) 
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where we have used Eq. (28), the relation la^°- = —1, and the fact that is a Kilhng 
vector. The result, of course, agrees with the standard one.|| 

It is also possible to show, using the expression for J"'' that the entropy in Eq. (93) 
is also equal to 

SNoether = ( T75 ) ^ab^cd.d'S, (96) 

where k is the surface gravity of the horizon and the [D — 2)-dimensional integral is on 
a spacelike bifurcation surface with eab denoting the bivector normal to the bifurcation 
surface. The variation is performed as if Rabcd and the metric are independent and the 
whole expression is evaluated on a solution to the equation of motion. A wide class 
of theories have been investigated using such a generalization in order to identify the 
thermodynamic variables relevant to the horizon. 



3. 6. The Lanczos-Lovelock models of gravity 

Among the class of theories described by the field equations 2Qab = Tab with Qab given 
by Eq. (81), one subset deserves special mention. These are the theories for which the 
Lagrangian satisfies the condition VaP"'^'^'^ = 0. (Since P"''^'^ has the symmetries of the 
curvature tensor, it follows that it will be divergence- free in all the indices.) In this case, 
Eq. (81) simplifies considerably and we get 

Gab = Pa^'^'^Rbcde — Tt^Qab = Tlab " -K^Qab] P"''"^'^ = — — (97) 

^ ^ O-ttabcd 

The crucial difference between Eq. (81) and Eq. (97) is the following. Since L does not 
depend on the derivatives of the curvature tensor, it contains at most second derivatives 
of the metric; therefore, P"-'"^'^ also contains only up to the second derivatives of the 
metric. It follows that Eq. (97) will not lead to derivatives of the metric higher than 
second order in the field equations. In contrast, Eq. (81) can contain up to fourth 
order in the derivatives of the metric. Though sometimes explored in the literature, 
field equations with derivatives higher than second order create several difficulties. (For 
example, when the equations are second order, the boundary condition in a variational 
principle is more easily defined than when higher order terms occur in the equations 
of motion; in such a case, the variation of the action functional requires very special 
procedures.) In view of this, there is a strong theoretical motivation to consider theories 
in which the Qab is of the form in Eq. (97). 

The Lagrangians which lead to the expression in Eq. (97) are, of course, quite 
special and are known as Lanczos-Lovelock Lagrangians. They can be expressed as a 

II There is, however, a subtlety that needs to be stressed regarding this derivation. In Einstein's theory, 
J° = 2i?^^^. Hence, for any static vacuum solution to Einstein's theory with being the Killing vector, 
the Noether current J° vanishes identically! A direct integration should therefore give zero entropy. 
This difficulty is circumvented by first obtaining J"'' and performing the integral over it on a single 
2-surface rather than integrating J° over a compact region in spacetime. The same situation arises in 
the calculation of Komar mass integrals for vacuum spacetimes. 
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sum of terms, each involving products of curvature tensors with the m— th term being 
a product of m curvature tensors. The general Lanczos-Lovelock Lagrangian has the 
form [66], 

^ 1 

^ - 2^ Cm-l-im) , ^(m) - "bib2...b2^ ^aiaa ^aa^-iaa™ ' 1^^^ 

m=l 

where the Cm are arbitrary constants and L(,m) is the m-th order Lanczos-Lovelock 
Lagrangian. The m = 1 term is proportional to S'^^R'^f oc R and leads to Einstein's 
theory. The m = 2 term gives rise to what is known as Gauss-Bonnet theory. Because 
of the determinant tensor, it is obvious that in any given dimension D we can only have 
K terms where D > 2K. It follows that, if D = 4, then only the K = 1,2 are non-zero. 
Of these, the Gauss-Bonnet term corresponding to -R' = 2 gives, on variation of the 
action, a vanishing bulk contribution. In dimensions D = 5 to 8, one can have both the 
Einstein-Hilbert term and the Gauss-Bonnet term etc. and so on for higher dimensions. 
It is conventional to take Ci = 1 so that the which gives Einstein gravity, reduces 
to (i?/167r). The normalizations m > 1 are somewhat ad-hoc for individual L(m) since 
the are unspecified at this stage. 

The Lanczos-Lovelock models possess black hole solutions and their thermodynamic 
properties have been investigated quite extensively. It can be shown, for example, that 
the entropy of a black hole horizon T-L in Lanczos-Lovelock models (determined using 
Eq. (96)) is given by 

S\n = ^Aiimcrn / rf^~^a;_Lv^^{m-i) , (99) 

m=l "^"^ 

where x_l denotes the transverse coordinates on "H, a is the determinant of the intrinsic 
metric on "H. It is interesting to observe that in these models, the entropy for the mth 
order theory is given by a surface integral involving the Lagrangian in the (m — l)th 
order theory. We will now indicate how this result arises. 

To do this we need to evaluate the Noether charge M corresponding to the 
current J", for a static metric with a bifurcation horizon and a Killing vector field 
= (1, 0); .^a = Qao- The location of the horizon is given by the vanishing of the norm 
^""^a = 900, of this Killing vector. The Noether charge is given by 



rf^-'xv^J°" (100) 

't,rH 

in which we have ignored the contributions arising from b when it ranges over the 
transverse directions. This is justifiable when transverse directions are compact or in 
the case of Rindler approximation when nothing changes along the transverse direction. 
In the radial direction, the integral picks out the contribution at r = rjy which is 
taken to be the location of the horizon. Using J"** = 2P"'"^"'Vc'Cd (see Eq. (86)) and 
= 9dai'' = 9 do we get: 

jor ^ 2pO^'='^V,ed = ^P'^^'^d.^d = 2P°^'=''9,to = ^P'^^'^ddg^o (101) 
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where we have used the symmetries of P"'^^'^ which are the same as those of the curvature 
tensor. So 



Jt,rH 

= 2m [ t/^-'xv^g^*'°9rf(7,o (102) 

Jt,rH 

where Q"''^'^'^ = (l/m)P"*^'^. Therefore the entropy is given by 

SNoether = I^U = 2/3m j ci^'^X yf^ Q^'^'^ddgcO (103) 

Jt,rH 

When the near- horizon geometry has a Rindler hmit, the r coordinate becomes the x 
coordinate and only goo = —k^x"^ contributes. Then this expression reduces to 

SNoether = 8nm [ rf^'^x^ (g°^) (104) 
Jh 

where a is the determinant of the metric in the transverse space. 

Let us consider this quantity Q^q for the m-th order Lanczos-Lovelock action, given 

by : 



QXO _ 1 1 rx0a3...a2m ( Tjbsb4 TDb2m~lb2m\ 



(105) 



where we have added a normahzation which gives Einstein's action for m = 1 and will 
define = l/lGvr for the m = case. We have also indicated that we are evaluating 
this expression in Rindler limit of the horizon, as in Eq. (104). The presence of both 
and X in each row of the alternating tensor forces all other indices to take the values 
2, 3, D - 1. In fact, we have CtX^ = ^PJttZ ^i^h A, = 2, 3, - 1 (the 
remaining combinations of Kronecker deltas on expanding out the alternating tensor are 
all zero since 5\ = Q = 5\ and so on). Hence reduces to 



_ \ (_}__}_ \ xA3^4...A2m ( T3B3B4, pB2„^-lB2„ 
° ~ 2 \ 167r2™"l ' "BsB4...B2m \,-^A3A4----^A2m-lA2„ 



(106) 



Therefore, in the e — )■ limit, recalling that |^= ^^Rcd Iw? "we find that Q^q is 
essentially the Lanczos-Lovelock Lagrangian of order (m — 1): 

Qlo = 2^{m-i) , (107) 

where we have restored the subscript giving the order of the Lagrangian. The entropy 
becomes 

4oLer = 4vrm [ rf^-2x^v^L(„_i) , (108) 
Jn 

This entropy in the m— th order Lanczos-Lovelock theory is an integral over the 
Lagrangian of (m — l)th order. For m = 1 (Einstein gravity), the L(o) is a constant 
giving an entropy proportional to transverse area; for m = 2 (Gauss-Bonnet gravity), 
the entropy is proportional to integral of R over transverse direction. 

While these results are satisfactory at a formal level, one must stress that the 
explicit value of the entropy depends on the nature of the theory decided by the 
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parameters Cj and may not have simple interpretation for certain range of parameters. 
For example, it is known that the on-shell entropy of the Lanczos-Lovelock models will 
not be positive definite for all range of parameters [67]. More seriously, the study of 
these models also raises several new conceptual conundrums for which it is difficult to 
find satisfactory answer. We shall now describe several of these issues. 

4. Thermodynamics of horizons: A deeper look 

4.I. The degrees of freedom associated with black hole entropy 

We have seen that there is a natural way of associating a temperature with any horizon 
including, for example, the Rindler horizon in fiat spacetime. But the arguments given in 
the last section leading to the association of entropy — in contrast to the association of 
a temperature — cannot be easily generalized from black hole horizon to other horizons. 
While there is general agreement that all horizons have a temperature, very few people 
[53, 68, 69, 70] have taken a firm stand as regards the question of associating an entropy 
with a horizon. To certain extent this ambivalence among researchers has led to most 
of the work being concentrated on analyzing black hole entropy (rather than horizon 
entropy) and we shall start our discussion with issues connected with black hole entropy. 

In the case of normal matter, entropy can be provided a statistical interpretation 
as the logarithm of the number of available microstates that are consistent with the 
macroscopic parameters which are held fixed. That is, S{E) is related to the degrees 
of freedom (or phase volume) g{E) by S{E) = In g{E). Maximization of the phase 
volume for systems which can exchange energy will then lead to equality of the quantity 
T{E) = {dS/dE)~^ for the systems. It is usual to identify this variable as the 
thermodynamic temperature. (This definition works even for self-gravitating systems 
in microcanonical ensemble; see eg., [71].) 

Assuming that the entropy of the black hole should have a similar interpretation, 
one is led to the conclusion that the density of states for a black hole of energy E = M 
should vary as 



Such a growth implies [69], among other things, that the Laplace transform of g{E) 
does not exist so that canonical partition function cannot be defined (without some 
regularization) . That brings us to the crucial question: What are the microscopic states 
which one should count to obtain the result in Eq. (109) ? That is, what are the degrees 
of freedom which lead to this entropy ? 

To begin with, the thermal radiation surrounding the black hole has an entropy 
which one can compute. It is fairly easy to see that this entropy will proportional to the 
horizon area but will diverge quadratically. Near the horizon the field becomes free and 
solutions are simple plane waves (see Section 2.5). It is the existence of such a continuum 
of wave modes which leads to infinite phase volume for the system. More formally, the 




(109) 
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number of modes n{E) for a scalar field with vanishing boundary conditions at two 
radii r = R and r = L is given by 



n{E) ~ — / ^ -il 

37r (1 - 2M/rf \ r 



3/2 



(110) 

in the WKB limit (see [72, 73]). This expression diverges as i? — )■ 2M showing that 
a scalar field propagating in a black hole spacetime has infinite phase volume. The 
corresponding entropy computed using the standard relations: 



S = p 



di3 



F = -rdE^^. (Ill) 



is quadratically divergent: S = {Ah/I'^) with / — )■ 0. The divergences described above 
occur around any infinite redshift surface and is a geometric (covariant) phenomenon. 

The same result can also be obtained from what is known as "entanglement entropy" 
arising from the quantum correlations which exist across the horizon. (For a review, 
see [74]). We saw in Section 3.2 that if the field configuration inside the horizon is 
traced over in the vacuum functional of the theory, then one obtains a density matrix 
p for the field configuration outside [and vice versa]. The entropy S = —Tr{p\n p) 
is usually called the entanglement entropy [75, 76, 77]. This is essentially the same 
as the previous calculation and, of course, S diverges quadratically on the horizon 
[78, 79]. In fact, much of this can be done without actually bringing in gravity into the 
picture; all that is required is a spherical region inside which the field configurations 
are traced out [80, 81]. Physically, however, it does not seem reasonable to integrate 
over all modes without any cut off in these calculations. By cutting off the modes at 
I Lp one can obtain the "correct" result but in the absence of a more fundamental 
argument for regularizing the modes, this result is not of much significance. The cut off 
can be introduced in a more sophisticated manner by changing the dispersion relation 
near Planck energy scales but again there are different prescriptions that are available 
[82, 83, 84, 85] and none of them are really convincing. 



4-.2. Black hole entropy in quantum gravity models 

There have also been attempts to compute black hole entropy in different models of 
quantum gravity [86]. In standard string theory this is done as follows: There are 
certain special states in string theory, called BPS states, that contain electric and 
magnetic charges which are equal to their mass. Classical supergravity has these states 
as classical solutions, among which are the extremal black holes with electric charge 
equal to the mass (in geometric units). These solutions can be expressed as a Reissner- 
Nordstrom metric with both the roots of Qqq = coinciding: obviously, the surface 
gravity at the horizon, proportional to gooirn) vanishes though the horizon has finite 
area. Therefore, these black holes have zero temperature but finite entropy. For certain 
compactification schemes in string theory (with d = 3,4,5 fiat directions), in the limit of 
G — 7- 0, there exist BPS states which have the same mass, charge and angular momentum 
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of an extremal black hole in d dimensions. One can explicitly count the number of such 
states in the appropriate limit and one finds that the result gives the density of states in 
Eq. (109) with correct numerical factors [87, 25, 88]. This is done in the weak coupling 
limit and a duality between strong coupling and weak coupling limits is used to argue 
that the same result will arise in the strong gravity regime. Further, if one perturbs the 
state slightly away from the BPS limit, to get a near extremal black hole and construct a 
thermal ensemble, one obtains the standard Hawking radiation from the corresponding 
near extremal black hole [88]. 

While these results are encouraging, there are several issues which are intriguing: 
First, the extremality (or near extremality) was used crucially in obtaining these results. 
We do not know how to address the entropy of a normal Schwarzschild black hole which 
is far away from the extremality condition. Second, in spite of significant effort, we do 
not still have a clear idea of how to handle the classical singularity or issues related to 
it. This is disappointing since one might have hoped that these problems are closely 
related. Finally, the result is very specific to black holes. One does not get any insight 
into the structure of other horizons, especially De Sitter horizon, which does not fit the 
string theory structure in a natural manner. 

The second approach in which some success related to black hole entropy is claimed, 
is in the loop quantum gravity (LQG). While string theory tries to incorporate all 
interactions in a unified manner, loop quantum gravity [89, 90] has the limited goal 
of providing a canonically quantized version of Einstein gravity. One key result which 
emerges from this program is a quantization law for the areas. The variables used in 
this approach are like a gauge field and the Wilson lines associated with them. The 
open Wilson lines carry a quantum number Jj with them and the area quantization law 
can be expressed in the form: Ah = SirG'y ^ Ji{Ji + 1) where Jj are spins defined on 
the links z of a spin network and 7 is free parameter called Barbero-Immirizi parameter. 
The Ji take half- integral values if the gauge group used in the theory is SU(2) and take 
integral values if the gauge group is SO (3). These quantum numbers, Ji, which live on 
the links that intersect a given area, become undetermined if the area refers to a horizon. 
Using this, one can count the number of microscopic configurations contributing to a 
given horizon area and estimate the entropy. One gets the correct numerical factor 
(only) if 7 = lnm/27r\/^ where m = 2 or m = 3 depending on whether the gauge group 
SU(2) or S0(3) is used in the theory [91, 92, 93, 94]. 

Again there are several unresolved issues. To begin with, it is not clear how exactly 
the black hole solution arises in this approach since it has been never easy to arrive at 
the low energy limit of gravity in LQG. Second, the answer depends on the Immirizi 
parameter 7 which needs to be adjusted to get the correct answer, after we know the 
correct answer from elsewhere. Even then, there is an ambiguity as to whether one 
should have SU(2) with 7 = In2/27ry2 or S0(3) with 7 = \n3/27iV2. The SU(2) was 
the preferred choice for a long time, based on its close association with fermions which 
one would like to incorporate in the theory. However, there has also been some occasional 
rethinking on this issue due to the following consideration: For a classical black hole. 
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one can define a class of solutions to wave equations called quasi normal modes [see 
e.g. ,[95, 96, 97, 98]]. These modes have discrete frequencies [99, 100, 101, 102] which 
are complex, given by 

n+(l/2) ln(3) , , 

The ln(3) in the above equation is not negotiable. If one chooses SO (3) as the gauge 
group, then one can connect up the frequency of quanta emitted by a black hole when the 
area changes by one quantum in LQG with the quasi normal mode frequency [103, 104]. 
It is not clear whether this is a coincidence or a result of some significance. If one assumes 
that this result is of some fundamental significance, the SO (3) gains preference. 

The short description given above shows that candidate models for quantum gravity 
are not yet developed well enough to provide a clear physical picture of horizon 
thermodynamics. (This is true even as regards numerous other approaches like e.g., 
those based on noncommutative geometry [105]). Given such a situation even in the 
well studied case of black hole horizon it is no suprise that we have virtually no quantum 
gravitational insight of other horizons. This fact gives additional impetus for studying 
the thermodynamic approach which we described in Section 1. There is, however, one 
central issue brought to the forefront by the quantum gravitational models of black hole 
entropy which we shall now discuss further. 



4-3. Black hole entropy: Bulk versus surface degrees of freedom 

Two obvious choices for the degrees of freedom contributing to the black hole entropy 
are those associated with the bulk volume inside the black hole (including those related 
to matter which collapses to form the black hole) or degrees of freedom associated with 
the horizon. One would have normally thought that the bulk degrees of freedom hidden 
by the horizon should scale as the volume V oc of the black hole. In that case, we 
would expect to get an entropy proportional to the volume rather than area. It is clear 
that, near a horizon, only a region of length Lp across the horizon contributes to the 
microstates so that in the expression {V/L^p), the relevant V is M'^Lp rather than M^. 
It is possible to interpret this as due to the entanglements of modes across the horizon 
over a length scale of Lp, which — in turn — induces a nonlocal coupling between the 
modes on the surface of the horizon. Such a field will have one-particle excitations, 
which have the same density of states as black hole [82, 83]. While this is suggestive 
of why we get the area scaling rather than volume scaling, a complete understanding is 
lacking. 

In fact, it is fairly easy to obtain an area scaling law for entropy if we assume that 
the degrees of freedom are on the horizon. Suppose we have any formalism of quantum 
gravity in which there is a minimum quantum for length or area, of the order of Lp. We 
have, for example, considerable evidence of very different nature to suggest Planck length 
acts as lower bound to the length scales that can be operationally defined and that no 
measurements can be ultra sharp at Planck scales [106]. Then, the horizon area Ah can 
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be divided into n = {Ah / CiLj:,) patches where ci is a numerical factor. If each patch has 
k degrees of freedom (due to the existence of a surface field, say), then the total number of 
microscopic states are fc" and the resulting entropy is S = nlnk = (41n fc/ci)(^ff/4Lp) 
which will give the standard result if we choose (41n /c/ci) = 1. The essential ingredients 
are only discreteness of the area and existence of certain degrees of freedom in each one 
of the patches. 

On the other hand, one cannot completely dismiss the degrees of freedom in the 
bulk as playing no role. Recall that the thermal density matrix for an eternal black 
hole can be obtained by explicitly integrating out the degrees of freedom on the left 
Rindler wedge (see Sec. 3.2). If integrating out certain degrees of freedom leads to a 
thermal density matrix, it makes sense to identify the same degrees of freedom as also 
contributing to the entropy. In which case, the bulk degrees of freedom has to play a 
role and still lead to an entropy that is proportional to the area. One possible way this 
can arise is through some kind of holographic relationship built in the theory of gravity 
so that the entropy of the bulk region V can be computed in terms of variables on its 
boundary dV. We shall see that the thermodynamic perspective of gravity, described 
later in this review, does lead to this possibility. 

4-4- Observer dependence of horizons and entropy 

We shall next address some conceptual issues brought about by the existence of horizons 
and entropy. We recall that the mathematical formulation leading to the association 
of temperature with any horizon is fairly universal and it does not distinguish between 
different horizons, like, for example, Rindler horizon in flat space or a Schwarzschild 
black hole event horizon or a de Sitter horizon. Assuming that temperature and 
entropy arise for fundamentally the same reason, it would be extremely unnatural not 
to associate entropy with all horizons. 

In this context, a distinction is sometimes made by arguing that, while the event 
horizon of a black hole can be given a purely geometrical definition, the Rindler horizon 
is observer dependent. This is irrelevant for the purpose of associating an entropy 
with the horizon. An observer plunging into a black hole will have access to different 
amount of information (and will attribute different thermodynamic properties to the 
black hole) compared to an observer who is remaining stationary outside the horizon. 
This is similar to what happens in the case of Rindler frame as well; an observer who 
stops accelerating or an inertial observer, will certainly have access to different regions 
of spacetime compared to the Rindler observer. In both the cases the physical effect of 
horizon in blocking information depends on the class of world lines one is considering. 
In that sense, all horizons are observer dependent. (There exists another example of a 
horizon — viz., the de Sitter horizon — the location of which depends on the observer 
but is considered "more real" than the Rindler horizon; those who object to assigning 
entropy to Rindler horizon are usually quite ambivalent about de Sitter or Schwarzschild- 
de Sitter [107] horizons!) It seems necessary to assign an (observer dependent) entropy 
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to all horizons. 

This feature, however, brings in a totally new layer of observer dependent 
thermodynamics into the theory which — though it need not come as a surprise — has 
to be tackled head-on. We know that an inertial observer will attribute zero temperature 
and zero entropy to the inertial vacuum. But a Rindler observer will attribute a finite 
temperature and non-zero (formally divergent) entropy to the same vacuum state. So 
entropy is indeed an observer dependent concept [108]. When one does quantum field 
theory in curved spacetime, it is not only that particles become an observer dependent 
notion but also the temperature and entropy. 

The observer in the Rindler wedge TZ will also perceive that the observables exhibit 
standard thermodynamic properties like entropy maximization, equipartition, thermal 
fluctuations etc., because the physics is governed by a thermal density matrix. But all 
these thermodynamical features arise because the Rindler observer attributes a density 
matrix to a pure quantum state after integrating out the unobservable modes. From 
this point of view, all these thermal effects are intrinsically quantum mechanical — 
which is somewhat different from the 'normal' thermal behaviour. Our results suggest 
that this distinction between quantum fluctuations and thermal fluctuations is artificial 
(like e.g., the distinction between energy and momentum of a particle in nonrelativistic 
mechanics) and should fade away in the correct description of spacetime, when one 
properly takes into account the fresh observer dependence induced by the existence of 
horizons. 

To see what all these imply in a concrete fashion, consider an excited state of a 
quantum field with energy 6E above the ground state in an inertial spacetime. When 
we integrate out the unobservable modes for the Rindler observer, we will get a density 
matrix pi for this state and the corresponding entropy will be 5*1 = — Tr (pi Inpi). The 
inertial vacuum state itself has the density matrix po and the entropy 5*0 = — Tr (po In Po) 
in the Rindler frame. The difference SS = Si — So is finite and represents the entropy 
attributed to this state by the Rindler observer. (This is finite though 5"! and 5*0 can be 
divergent.) In the limit of k — )■ cxd, which would correspond to a Rindler observer who 
is very close to the horizon, we can actually compute it and show that 

SS = (3SE = —5E (113) 

K 

To prove this, note that if we write pi = Po + Sp, then in the limit of k — )■ oo we can 
concentrate on states for which 5p/po ^ 1- Then we have 

-6S = Tr (pi In pi) - Tr (po In po) ^ Tr {6p In po) 

= Tr {6p{-I3Hr)) = -/3Tr ((pi - po)Hr) = -/SSE (114) 

where we have used the facts Tr (5p ~ and po = exp{—l3Hfj) where Hfj is the 
Hamiltonian for the system in the Rindler frame. The last line defines the 6E in terms 
of the difference in the expectation values of the Hamiltonian in the two states. This 
is the amount of entropy a Rindler observer would claim to be lost when the matter 
disappears into the horizon. (This result can be explicitly proved for, say, one particle 
excited states of the field [109].) 
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The above result is true in spite of the fact that, formally, matter takes an infinite 
amount of coordinate time to cross the horizon as far as the outside observer is concerned. 
This is essentially because we know that quantum gravitational effects will smear the 
location of the horizon by 0{Lp) effects [106]. So one cannot really talk about the 
location of the event horizon ignoring fluctuations of this order. From the operational 
point of view, we only need to consider matter reaching within few Planck lengths of the 
horizon to talk about entropy loss. In fact, physical processes very close to the horizon 
must play an important role in order to provide a complete picture of the issues we are 
discussing. There is already some evidence [82, 83] that the infinite redshift induced 
by the horizon plays a crucial role in this though a mathematically rigorous model is 
lacking. 

One might have naively thought that the expression for entropy of matter crossing 
the horizon should consist of its energy 5E and its own temperature Tmattcr rather than 
the horizon temperature. But the correct expression is 5S = SE/T^orizon] the horizon 
acts as a system with some internal degrees of freedom and temperature Thorizon as far 
as Rindler observer is concerned so that when one adds an energy 6E to it, the entropy 

change is 6S = ( / Thorizon )• 

All these are not new features but only the consequence of result that a Rindler 
observer attributes a non-zero temperature to inertial vacuum. This temperature 
influences every other thermodynamic variable. Obviously, a Rindler observer (or an 
observer outside a black hole horizon) will attribute all kinds of entropy changes to 
the horizons she perceives while an inertial observer (or an observer falling through the 
Schwarzschild horizon) will see none of these phenomena. This requires us to accept 
the fact that many of the thermodynamic phenomena needs to be now thought of as 
specifically observer dependent. 

5. Action functionals for gravity and horizon thermodynamics 

While deriving the expression for the temperature associated with a horizon we stressed 
the fact that it was completely independent of the dynamical equations satisfied by 
the metric. For example, we never needed to use Einstein's equations in the case of 
Schwarzschild black hole, say, to obtain the expression for temperature. The situation 
as regards Hawking evaporation is also similar; given a specific form of the metric, 
one could obtain Hawking evaporation without demanding that the metric should be a 
solution to any specific field equation. 

It therefore comes as a surprise that there is a deep and curious connection between 
the field equations of gravity and the horizon thermodynamics. We shall first provide a 
simple illustration of this result and will then consider a more general approach. 
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5.1. An unexplained connection between horizon thermodynamics and gravitational 
dynamics 

To illustrate this result in the simplest context [38], let us consider a static, spherically 
symmetric horizon, in a spacetime described by a metric: 

ds^ = -f{r)^dt' + f~\r)dr^ + r^d^^. (115) 

(These results can be easily generalized to the case with qqq ^ —g^^ by changing r^dVl^ 
by R^{r)dVt^ where R{r) is an arbitrary function. We will not bother to do this.) Let 
the location of the horizon be given by the simple zero of the function /(r), say at r = a. 
The Taylor series expansion of /(r) near the horizon /(r) ^ /'(a)(r — a) shows that the 
metric reduces to the Rindler metric near the horizon in the r — t plane with the surface 
gravity k = (c^/2)/'(a). Then, an analytic continuation to imaginary time allows us to 
identify the temperature associated with the horizon to be 

hcf\a) 

knT = (116) 

where we have introduced the normal units. The association of temperature in Eq. (116) 
with the metric in Eq. (115) only requires the conditions /(a) = and /'(a) 7^ 0. The 
discussion so far did not assume anything about the dynamics of gravity or Einstein's 
field equations. 

We shall now take the next step and write down the Einstein equation for the metric 
in Eq. (115), which is given by (1 — /) — rf'{r) = — (87rG/c^)Pr^ where P = is the 
radial pressure. When evaluated on the horizon r = a we get the result: 



inPa^ (117) 



G 

If we now consider two solutions to the Einstein's equations differing infinitesimally 

in the parameters such that horizons occur at two different radii a and a + da, then 

multiplying the Eq. (117) by da, we get: 
4 4 

^f'{a)ada - ^da = P{4na^da) (118) 
2G 2G 

The right hand side is just PdV where V = (47r/3)a^ is what is called the areal volume 
which is the relevant quantity when we consider the action of pressure on a surface area. 
In the first term, we note that /'(a) is proportional to horizon temperature in Eq. (116). 
Rearranging this term slightly and introducing a h factor by hand into an otherwise 
classical equation to bring in the horizon temperature, we can rewrite Eq. (118) as 





2 ^ Pd\—a'] (119) 



dS PdV 
The labels below the equation indicate a natural — and unique — interpretation for 
each of the terms and the whole equation now becomes TdS = dE + PdV allowing us 
to read off the expressions for entropy and energy: 

T IAh ^ { Ah^ ^'^ 



S = ^(ATTa') = E = -—a=-{ — ] (120) 
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where is the horizon area and Lp = Gh/c'. The result shows that Einstein's 
equations can be re-interpreted as a thermodynamic identity for a virtual displacement 
of the horizon by an amount da. 

The uniqueness of the factor P{4:Tia^)da, where Atto^ is the proper area of a surface 
of radius a in spherically symmetric spacetimes, implies that we cannot carry out the 
same exercise by multiplying Eq. (117) by some other arbitrary factor F{a)da instead of 
just da in a natural fashion. This, in turn, uniquely fixes both dE and the combination 
TdS. The product TdS is classical and is independent of h and hence we can determine 
T and S only within a multiplicative factor. The only place we introduced h by hand 
is in using the Euclidean extension of the metric to fix the form of T and thus 5*. The 
fact that T (X h and S (xl/his analogous to the situation in classical thermodynamics 
in contrast with statistical mechanics. The TdS in thermodynamics is independent of 
Boltzmann's constant while statistical mechanics will lead to S (x ks and T oc l/ks- 

It must be stressed that this result is quite different from the conventional first 
law of black hole dynamics, a simple version of which was mentioned in Eq. (2). The 
difference is easily seen, for example, in the case of Reissner-Nordstrom black hole 
for which = P is non-zero due to the presence of nonzero electromagnetic energy- 
momentum tensor in the right hand side of Einstein's equations. If a chargeless particle 
of mass dM is dropped into a Reissner-Nordstrom black hole, then the standard first 
law of black hole thermodynamics will give TdS = dM. But in Eq. (119), the energy 
term, defined as E = a/2, changes by dE = {da/2) = (1/2) [a/(a - M)]dM ^ dM. It 
is easy to see, however, that for the Reissner-Nordstrom black hole, the combination 
dE + PdV is precisely equal to dM making sure TdS = dM. So we need the PdV term 
to get TdS = dM from Eq. (119) when a chargeless particle is dropped into a Reissner- 
Nordstrom black hole. More generally, if da arises due to changes dM and dQ, it is 
easy to show that Eq. (119) gives TdS = dM — {Q/a)dQ where the second term arises 
from the electrostatic contribution. This ensures that Eq. (119) is perfectly consistent 
with the standard first law of black hole dynamics in those contexts in which both are 
applicable but dE ^ dM in general. It may also be noted that the way Eq. (119) was 
derived is completely local and quite different from the way one obtains first law of black 
hole thermodynamics. 

It is quite surprising that the Einstein's field equations evaluated on the horizon 
reduces to a thermodynamic identity. More sharply stated, we have no explanation as to 
why an equation like Eq. (119) should hold in classical gravity, if we take the conventional 
route. This strongly suggests that the association of entropy and temperature with a 
horizon is quite fundamental and is actually connected with the dynamics (encoded in 
Einstein's equations) of the gravitational field. The fact that quantum field theory in 
a spacetime with horizon exhibits thermal behaviour should then be thought of as a 
consequence of a more fundamental principle. 
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5.2. Gravitational field equations as a thermodynamic identity on the horizon 

If this conjecture is correct, the equahty — between field equations on the horizon 
and the thermodynamic identity — should have a more general validity. This 
has now been demonstrated for an impressively wide class of models like the cases 
of stationary axisymmetric horizons and evolving spherically symmetric horizons in 
Einstein gravity [110], static spherically symmetric horizons [111] and dynamical 
apparent horizons [112] in Lanczos-Lovelock gravity, generic, static horizon in Lanczos- 
Lovelock gravity [113], three dimensional BTZ black hole horizons [114], FRW 
cosmological models in various gravity theories [115, 116, 117, 118, 119, 120, 121, 122, 
123] and even in [124] the case Horava-Lifshitz gravity. 

We shall describe briefly how this is achieved in the case of an arbitrary, static, 
horizon in Einstein's theory and in Lanczos-Lovelock models (for more details, see ref. 
[113]). Consider a static spacetime with the metric 

ds^ = -NMt^ + dn^ + aABdy^dy^ (121) 

where aAsin^y^) is the transverse metric, and the Killing horizon, generated by the 
timelike Killing vector field ^ = dt, is approached as A^^ — )■ 0. Near the horizon, 
~ + 0{n^) where k is the surface gravity [125]. The t = constant part of the 
metric is written by employing Gaussian normal coordinates for the spatial part of 
the metric spanned by (n, y^) with n being the normal distance to the horizon. By 
manipulating the Einstein's equations evaluated on the horizon, one can prove [113] the 
following relation: 

= r|V^M (122) 

where A is the affine parameter along the outgoing null geodesies and the R\\ is the Ricci 
scalar of the on- horizon transverse metric, [ct//]^^. The Einstein tensor components are 
evaluated in an orthonormal tetrad appropriate for a timelike observer moving along 
the orbit of the Killing vector field generating the Killing horizon. This is denoted by a 
hat on the indices; for example, ^ = {—gtt)"^^'^ dt etc., and — G| = G^= G{^,^). We 

have used \h = \h in the second equality in Eq. (122) and Einstein's equation in 

the third one. (The fact that Cl = on the horizon is crucial; more details regarding 
this symmetry can be found in Ref. [125].) 

Multiplying Eq. (122) by d^y, and integrating over the horizon 2-surface, we obtain 

a r n _ .o 1 r r 1 „ ^ .o i a 
d\ 



J 


f\y^d'y 


6X- 








H 





[ P^^d^ySX (123) 

JH 



where we have identified T = h/2t[ as the horizon temperature, and used the 
interpretation of as the normal pressure. Pi, on the horizon. We can therefore 
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interpret 

F = [ P^V^d'y (124) 
Jh 

as the the average normal force over the horizon "surface" and F 5X as the (virtual) 
work done in displacing the horizon by an affine distance 6X. Equation (123) can now 
be written as 

T5xS -5xE = F SX (125) 

where 

S=-^j^<l'y (126) 

is (a priori) just a function of A; in particular, the derivative of S with respect to A is 
well-defined and finite on the horizon. We only need the expression for S very close to 
the horizon. The value of S" at A = Xh-, 

S{X = Xh) = \I v^d^y (127) 

is equal to the Bekenstein-Hawking entropy of the horizon. We also identify the energy 
E associated with the horizon as 

where x is the Euler characteristic of a 2-dimensional compact manifold A4.2 which in 
this case would be the horizon 2-surface, given by 

X {M2) = ^ [ R d[vol] (129) 
4vr Jm2 

(If the manifold has a boundary, then the expression for Euler characteristic will have 
additional boundary terms.) Thus Einstein's equations evaluated on the horizon can be 
expressed as a thermodynamic identity and — as a bonus — we get a geometric definition 
of energy. Our particular identification of E is fixed by the choice of the affine parameter 
along the outgoing null geodesies. In particular, this brings out the significance of the 
radial coordinate r in spherically symmetric and stationary spacetimes; in either case, 
r is the affine parameter along the outgoing null geodesies. 

To connect up with the previous discussion, let us consider again the spherically 
symmetric case with a compact horizon, in which A = r and x = 2. We obtain E = rnfi, 
with rn being the horizon radius, which matches with the standard expression for 
quasilocal energy for such spacetimes obtained previously. In general, for a compact, 
simply connected horizon 2-surface, X = 2 (since any such manifold is homomorphic to 
a 2-sphere), and we have, E = Xh /2. Therefore, for spherically symmetric black holes, 
since Pi = Pr is independent of the transverse coordinates {9,(f)), we obtain 

T6S -6E = Pr6V (130) 

with 6S = 27rrHSrH, SE = drn/^, Pr = T^i^rn), SV = Anrj^ dru and T is the horizon 
temperature. We therefore recover the result in Eq. (119). 
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Exactly similar structure emerges for the near-horizon field equations of Lanczos- 
Lovelock gravity as well. In this case, the analysis proceeds along identical lines 
though the algebra is more complicated. To begin with, using the field equation 
= {\/2)T^ on the horizon for the mth order Lanczos-Lovelock theory in D dimensions 
and manipulating the expressions, one can obtain [113] the relation: 

- L(^-')v^ 6\ + 0[i\ - XhY^' 6X] (131) 

where T is the horizon temperature, Lm^ is the Lanczos-Lovelock Lagrangian for the 
mth order Lanczos-Lovelock theory in D dimensions and 

_ rBAi..._Bm_i {D-2)pCiDi {D-2) pCm^iDm-i ( ^ "iO^ 

- OACi...Dm-i ^ AiBi ■ ■ ■ A„_iB™_i l-^-^ZJ 

where the upper case Latin indices A,B,.... etc. run over the transverse coordinates. 
We can now prove that the factor multiplying T in Eq. (131) is directly related to the 
variation of the following quantity, with the variation being evaluated at A = Xh'- 

S = Anm [ dS L^^S^^^ (133) 



To do this, we use the fact that the variation of S in Eq. (133) must give equations 
of motion for the (m — 1)"^ order Lanczos-Lovelock term in {D — 2) dimensions. (The 
variation would also produce surface terms, which would not contribute when evaluated 
at A = Xh because the horizon is a compact surface with no boundary.) We therefore 
have: 

5^S = -47rm [ dS £'^^ Sxctbc (134) 
Jn 

where the variation is evaluated on A = Xh- Noting that the Lagrangian is 



and 



r{D-2) _ 1 1 (D-2)pC m — 1 

— -^Q^ 2{m-l) C1D1...D 

m — 1 -^m — 1 — 1 



1 1 1 

r>i i^iaibi...ambm -pcidi 

TDCmdm 

aifei 



2 IGvr 2™ icidi...Cmdm aifei ambm 

——ST^f'-T^rR'th ■■■RYh --5]Lm (135) 

IGTT 2™ dl...Cmdm aiOi ambm 2 ^ 



we see that 



S'^ = -- — -^S^ (136) 

^ 2 167r2('"-i) ^ ^ ' 



Therefore, we obtain 



i m 

0\b 



^ /" dS 6xaBC (137) 



^ 2(™- 

which is precisely the integral of the factor multiplying T in Eq. (131). As mentioned 
earlier, S defined in Eq. (133) is a function of A, and its derivative with respect to A is 
well defined and finite on the horizon. The expression for 5*, evaluated at A = Xh-, 



S{X = Xh) = 47rm / dS L[^_i ^ (138) 
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can be interpreted as the entropy of the horizon and it matches with the standard result 
in Eq. (108) obtained by other methods. Multiplying Eq. (131) by S^~'^^y, integrating 
over the horizon surface, and taking the n — j- limit, we now see that it can be written 

as 

TSxS- I ds s\= / ds ri s\= I dE r| s\ 



' n 



H JH ^ JH 

dS 6X (139) 

H 

where we have used the field equations = (1/2)T| in the first equality, and the 

relation \h = \h the second equality. This equation now has the desired form 
of the first law of thermodynamics, when we identify: (i) the S, defined by Eq. (138), as 
the entropy of horizons in Lanczos-Lovelock gravity; indeed, exactly the same expression 
for entropy has been obtained in the literature using the Wald entropy (e.g, ref. [126]) 
and (ii) the second term on the left hand side of Eq. (139) as S\E] this leads to the 
definition of E to be 

E = J\\ J dS (140) 

where the A — )■ Xh limit must be taken after the integral is done. (Therefore, we need 
to know the detailed form of L function of A to calculate this explicitly). For 

D = 2(m + 1), the integral over "H in Eq. (140) is related to the Euler characteristic 
of the horizon and we get E oc A^. For m = 1, D = 4, this reduces to the expression 
obtained earlier in the case of Einstein gravity. 

One can also determine the scaling of E for spherically symmetric spacetimes for 
general Lanczos-Lovelock lagrangians, with horizon at r = rj^, and A = r. In this case, 
L^f"^^ oc (l/A^)™ and ^/^ oc X^'^. The integrand in Eq. (140) scales as \iD~2)~2m 
giving E oc A^ As mentioned above, for D = 2(m + l), E oc A^- In fact, 

in the case of spherically symmetric spacetimes in Lanczos-Lovelock theory, the above 
expression can be shown to be exactly equivalent to the one derived by others (see [111], 
and also Ref. [14] therein). No general expression for energy in Lanczos-Lovelock theory 
exists in the literature, and Eq. (140) could be thought of as first such definition which 
appears to be reasonable from physical point of view. 

We thus conclude that for generic static spacetimes in Lanczos-Lovelock gravity, 
the field equations can be written as a thermodynamic identity: 

T6xS - 6xE = F6X (141) 

thereby showing that the thermodynamic relations hold for much more general cases 
than Einstein's theory. 



5. 3. Structure of gravitational action functionals 

Since the field equations of a theory can be derived by varying the dynamical variables 
in a suitably defined action functional, it makes sense to examine the nature of 
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the action functionals in order to learn more about the curious connection between 
horizon thermodynamics and gravitational dynamics discussed above. We shall now 
describe several peculiar features of action functionals in Einstein's theory as well as in 
more general Lanczos-Lovelock theories of gravity. Once again, these results have no 
explanation in the conventional approach. 

Consider a theory based on a gravitational Lagrangian Lg which can be expressed 
in the form 

IQnLg ^ Q:^'R\,, = QX (142) 

where the tensor Qj"^'^ has (i) all the symmetries of the curvature tensor and (ii) has zero 
divergence on all the indices, VaQ"'"^'^ = etc. With some simple algebraic manipulation, 
we can express [41] any such action in the form 

= V^-^quad + -^sur (143) 

where we have separated out the expression into one term [Lquad] which is quadratic in 
(and hence quadratic in the first derivatives of the metric) and another term [Lguv] 
which is a total divergence that could lead to a surface term on integration. Thus, to 
begin with, we have the result that a wide class of theories determined by the tensor 
Qabcd ]^g^g action which will naturally lead to a surface term. We shall next show 
that: (i) there is a holographic relationship between the surface and bulk term and that 
(ii) the surface term leads to the entropy of the horizon. 

5.4- Holographic nature of surface term in Einstein- Hilbert action 

To set the stage, let us begin with Einstein's theory which has the Lagrangian Lg = 
(167r)~^i?. This Lagrangian can be expressed in the form in Eq. (142) with 

Qa^ = li^^g'' - 5^9'% Qab = = - 5t5l) (144) 

In this particular case, the bulk term involving the quadratic part is given by 

I^quad = 2Q^"^'T'^,,Tl = g^" (T]Ji - T^^l^ (145) 
and the surface term arises from 

Lsur = 2d, [^Qj>^ri] = 2d, [V^Q^ig'Tl] ^ d, [y^V^] (146) 
where we have defined a four component object (which is not a four-vector) by: 

^ {g^'Vr, - /'^r™ ,) = ((7-/^- - g^'g")d,gaj = --^d,{gg'^) (147) 

By direct computation one can verify that the bulk and the surface term are related by 
the equation 

y/^Lsur = -da \9ij ^^(^^ ^Y^ (^'^^) 



didagij) 
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(We call this relation 'holographic' for want of better terminology; it has no connection 
with the term holography used in string theory.) This structure has a simple physical 
meaning that can be understood as follows [39, 40]. Given a Lagrangian Lq{q,q) in 
classical mechanics, say, one can obtain the standard Euler-Lagrange equations by 
varying q in the action functional with the condition 5g = at the boundary. Consider 
now a different Lagrangian defined as 



(149) 



This Lagrangian, unlike Lq contains q. If we vary the action resulting from Lp{q,q,q) 
but — instead of demanding 6q = at the boundary — demand that 6p = at 
the boundary where p{q,q) = [dLg/dq) is the momentum, then we will get the same 
equations of motion as the one obtained from Lg. That is, even though Lp contains the 
second derivatives of q, it leads to second order differential equations for q (rather than 
third order) if we fix p at the boundary. Lagrangians involving second derivatives of 
dynamical variables but in a specific combination through the second term in Lp{q, q, q) 
are quite special. This idea generalizes trivially to field theory and we see, comparing 
Eq. (148) with Eq. (149), that Einstein's theory has this special structure. It is this 
holographic relationship which allows the surface terms to contain information about 
the bulk. 

5.5. Horizon entropy and the surface term in Einstein- Hilbert action 

It is also easy to show that the surface term actually leads to the horizon entropy in the 
case of Einstein's theory. To do this, we will work with the Euclidean extension of the 
action in which the horizon is mapped to the origin. Close to the origin, in Rindler like 
coordinates, we have the metric of the form 

ds"^ = n'^^'^dtl + + dyi\ (150) 

To evaluate the surface integral arising from Lgur in Eq. (146) on the horizon, we shall 
compute it on the surface ^ = e around the origin in the ^ — t_E plane and then take 
the limit of e — )■ 0. So we need to integrate \/hncV'^ where is defined by Eq. (146), 
= i%ey/a with a being the determinant of the metric in the transverse coordinates 
and n* = 5^ is the normal. In the integral, the range of tE-, being an angular coordinate, 
is (0, /3 = 2tx/k). Using = — 2/e, we find that the surface contribution to the action 
is 

167rAur = J d^xVhncV = J dts J d'^:s.±{Key/a) ^-'^ 

= -47rA, (15V 



Therefore 



Aur = -jA^ (152) 
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where A± is the transverse area. This resuh shows that the surface term in the action 
has a direct thermodynamic meaning as the horizon entropy [127]. The sign flips if we 
change the sign of ric and hence is not of real significance. (But, with our choice for Uc, 
the sign can be explained by the fact that the probability for a configuration is related 
to Euclidean action hj V oc exp(— ^sur) while V cc exp(5') where S is the entropy; 
hence S = — ^sur-) As we said before, this result has no explanation in the conventional 
approach in which the field equations know nothing about the surface term. 

There is another curious aspect related to the surface term in Einstein-Hilbert 
action which is worth mentioning. In standard quantum field theory, the kinetic term 
for the field (p will be quadratic in the derivatives of the field variable (((90)^) which 
will be integrated over the four volume to obtain the action. In natural units, action 
is dimensionless and hence all fields will have the dimension of inverse length. In the 
case of gravitational field, one might like to associate a second rank symmetric tensor 
field, Hab, to describe the graviton. In that case, the metric gab will be interpreted as 
Qab = Vab + ^-f^ab whcrc / is a coustaut with dimensions of length. [In normal units, 
P oc (Gh/c^).] Let us consider what happens to the Einstein-Hilbert action when we 
use this expansion and retain terms up to the lowest non-vanishing order in the bulk and 
surface terms. Expanding Lquad and Lsur in Taylor series in / and choosing P = IQtcG, 
where G is the Newtonian gravitational constant re-introduced for the sake of clarity, 
we find that the action functional becomes 

A = — ^ / d^Xy/^R = ^quad + Am- (153) 



where 



and 



with 



Aw = \j M'^'''^\r^^^)daHb,d,H,k + 0{l) (154) 
Aur = Jijd^^ dadb[H^' - v'^Hi] + 0(1) (155) 



This Auad matches exactly with the action for the spin-2 field known as Fierz-Pauli 
action (see e.g. Ref. [128]). However, the surface term — which is usually ignored in 
the discussions — is non-analytic in the coupling constant. Hence one cannot provide 
an interpretation of black hole entropy (which, as we have seen, can be obtained from 
the surface term in the action) in the linear, weak coupling limit of gravity. The integral 
we evaluated in the Euclidean sector around the origin to obtain the result in Eq. (152) 
cannot even be defined usefully in the weak field limit because we used the fact that goo 
vanishes at the origin. When we take goo = Voo + ^oo and treat hoo as a perturbation, it 
is obviously not possible to make goo vanish. 

In fact, the non- analytic behaviour of Aur on / can be obtained from fairly simple 
considerations related to the algebraic structure of the curvature scalar. In terms of 
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a spin-2 field, the final metric arises as gab = f]ab + ^ Hat where / oc ^/G has the 
dimension of length and Hab has the correct dimension of (length)"^ in natural units with 
h = c = 1. Since the scalar curvature has the structure R ~ (dg)"^ + d'^g, substitution 
of gab = Vab + I Hab gives to the lowest order: 

Leh oc ^i? ~ (dHf + jd'H (157) 

Thus the full Einstein-Hilbert Lagrangian is non-analytic in / because of the surface 
term. 



5. 6. Holographic structure of action functional in Lanczos-Lovelock gravity 

We shall now consider the generalization of these results to Lanczos-Lovelock theories 
with an action described by a tensor Q"-^'^'^. In this case, Q"'^'^'^ depends on the metric 
as well as curvature tensor but not on the derivatives of the curvature tensor. The 
Lagrangian for the m— th order Lanczos-Lovelock theory in D— dimension is given by 

T, ^ _ rl357...2fc-lp24p68 p2A:-2 2fc . t. _ n f^KS^ 

where k = 2m is an even number. The is clearly a homogeneous function of degree 
m in the curvature tensor i?"^ so that it can also be expressed in the form: 



^(m) 
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dL 



(m.) 
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PI bed T->a 



bed 



m 

jabed ^/^abed 



bed' 



(159) 



where P/^'^ = ((9L(^)/9i?'^ft^J so that P""^'' = mQ" 

The canonical momentum conjugate to the metric has to be defined more carefully 
in this case because Lbuik will contain second derivatives of the metric (unlike in the 
case of Einstein's theory). Once this technical problem is taken care of, one can show by 
direct computation that the surface and bulk terms obey a holographic relation given 
by 



[{D/2) - m]L(-) = -d. 



9ab 



"^bulk 



dig. 



dL 



5{digab) ' "^'dididjgab) 
dK[(t>,di(t>,.. 



(m) 
bulk 



where the Euler derivative is defined as 



da 



+ 



(160) 



(161) 



d{da(t>) 

The proof involves straight forward combinatorics [41]. This shows that a wide class of 
gravitational theories which have the surface term in the action functional exhibits the 
holographic relationship between surface and bulk terms. 



5. 7. Horizon entropy from the surface term in the Lanczos-Lovelock action functional 

To generalize the result that the surface term in the action functional gives the entropy 
of the horizon, we need to compare the surface term with the entropy in the Lanczos- 
Lovelock theories which was obtained earlier in Eq. (96) by using the Noether charge. 
We saw that, the Noether charge approach leads to the expression (see Eq. (103)): 

SNoether = / d^^^^X ^ Q"''° d^g^o', P = — (162) 

Jt,rH ^ 
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We will now show that the same result can be obtained by evaluating the surface term 
in the action on the horizon. For the Lanczos-Lovelock models, the bulk and the surface 
Lagrangians are given by 

i^quad = 2Qj"^'T'^,,Tl, Ls.,. = 2d, [V^Q^^Xd] (163) 
In the stationary case, the contribution of surface term on the horizon is given by 

= 2 I dt I d''-'x./^Q'''"^%gad (164) 

Once again, taking the integration over t to be in the range (0, /3) and ignoring transverse 
directions, we get 

Ssur = 2f3 [ rf^-^x Q^'^'d^gao (165) 



Comparing with Eq. (162), we find that 

SNoether = raSgur (166) 

The overall proportionality factor m has a simple physical meaning. Eq. (160) tells us 
that the quantity mLsur, rather than Lgur, which has the ^^d{qp)" structure and it is 
this particular combination which plays the role of entropy, as to be expected. 

It must be stressed that the above results defy understanding in the conventional 
approach. To begin with, it is not clear why the simplest generally covariant action in 
general relativity (and in Lanczos-Lovelock models) contains a total divergence term 
which leads to a surface term. Further, in the conventional approach, we ignore the 
surface term completely (or cancel it with a counter-term) and obtain the field equation 
from the bulk term in the action. Any solution to the field equation obtained by this 
procedure is logically independent of the nature of the surface term. But when the 
surface term (which was ignored) is evaluated at the horizon that arises in any given 
solution, it does correctly give the entropy of the horizon! This is possibly because of the 
specific relationship between the surface term and the bulk term given by Eq. (148) and 
Eq. (160). But these relations are again totally unexplained feature in the conventional 
approach to gravitational dynamics. Given that the surface term has the thermodynamic 
interpretation as the entropy of horizons, and is related holographically to the bulk term, 
we are again led to an indirect connection between spacetime djTiamics and horizon 
thermodynamics. 

5.8. Gravitational action as the free energy of spacetime 

There is one more aspect of the gravitational action functional which adds strength to 
the thermodynamic interpretation. We can show that, in any static spacetime with 
a bifurcation horizon (so that the metric is periodic with a period /3 = 27r/K in the 
Euclidean sector), the action functional for gravity can be interpreted as the free energy 
of spacetime. 
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Since the spacetime is static, there exists a timehke Kilhng vector field with 
components = (1, 0) in the natural coordinate system which exhibits the static nature 
of the spacetime. The conserved Noether current for the displacement x"" x"" + C,"" is 
given by Eq. (89) . We will work in the Euclidean sector and integrate this expression, 
taken on-shell with 2Qab = Tab, over a constant-t hypersurface with the measure 
dTja = 5^N\/hd^~^x where Qqq = N"^ and h is the determinant of the spatial metric. 
Multiplying by the period /3 of the imaginary time, we get 



/3 



J rdEa = (3 j T^edT^a + 13 j LCdT^a 

= j {l3N)T^UaU^^d^-^x + dtE j L^d^'^x (167) 

where we have introduced the four velocity u"" = ^""/N = N~^6q of observers moving 
along the orbits of and the relation dJ^a = Ua\/hd^~^x. The term in Eq. (106) 
involving the Lagrangian gives the Euclidean action for the theory. In the term involving 
Tab we note that /3iV = [3\oc corresponds to the correct redshifted local temperature. 
Hence we can define the (thermally averaged) energy E as 

j {PN)T^UaU^Vhd^~^x = J PiocTbUaU^Vhd'^^^x = l3E (168) 

We thus get 

A = (3 j rdT.a-l3E (169) 

We have, however, seen earlier that the first term involving the Noether charge gives 
the horizon entropy, which continues to hold true in the Euclidean sector. Therefore, 
we find that 

A = S-PE = -(3F (170) 

where F is the free energy. (Usually, one defines the Euclidean action with an extra 
minus sign a.s A = —Ae in which case the Euclidean action can be interpreted directly 
as the free energy.) The motivation for the definition of E in Eq. (168) now becomes 
clear since the entropy is related to the spatial integral of the energy density p(x) with a 
weightage factor /3(x) when the temperature varies in space. One can also obtain from 
Eq. (167) the relation: 

S = (3 [ rdT.a= I ^/^d^x{p + L) (171) 



where p = TabU°'u^ ■ This equation gives the entropy in terms of matter energy density 
and the Euclidean action. Alternatively, if we assume that the Euclidean action can be 
interpreted as the free energy, then these relations provide an alternate justification for 
interpreting the Noether charge as the entropy. (Similar results have been obtained in 
Ref. [129] by a more complicated procedure and with a different motivation.) 

There is another result which one can obtain from the expression for the Noether 
current. Taking the component of Eq. (89) and writing J° = Vfc we obtain 

^ =94v^J°")-2^° (172) 



-9 
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Only spatial derivatives contribute in the first term on the right hand side when the 
spacetime is static. This relation shows that the action obtained by integrating Ly/—g 
will generically have a surface term related to J"^ (In Einstein gravity Eq. (172) will 
read as L = 2Rq — 2Gq; our result generalizes the fact that Rq can be expressed as 
a total divergence in static spacetimes [130].) This again illustrates, in a very general 
manner, why the surface terms in the action functional lead to horizon entropy. 

Eq. (172) can be expressed more formally by introducing a vector = V^t = 
dkt = (1,0) which is the unnormalised normal to t = constant hypersurfaces. The unit 
normal is Ik = Nl^ = —Uk in the Lorentzian sector. Contracting with the expression 
for Noether charge, written for a Killing vector in the form 

2gl;^'' + Li^ = .]'' = VaJ'"' (173) 
and using Ik^ aJ'^"' = Va(/fcJ*^°), we get 

L = VaihJ'n - 2gU\ = -^d^ {V^ /fe J^") - 2g^eik (174) 

which is identical to Eq. (172) but is sometimes more convenient for manipulation. 

Finally, we show how several of the results obtained in ref.[69] giving a 
thermodynamic interpretation for Einstein-Hilbert action, can be generalized for 
Lanczos-Lovelock gravity. Using the Lanczos-Lovelock equations of motion, we can 
easily show that, for the m-th order Lanczos-Lovelock model in D-dimensions, 

T = -{D-2m)L; 20^ + L6^, = T,^ - -^^^—^6^,T (175) 
Defining as before p = TabU°'u^ and 

we have the relations 

(e - p) = -L; VaihJ'n = -e (177) 

Integrating over spacetime, treating the time integral as multiplication by /3, and defining 
the entropy in terms of the Noether charge, we get 

S - I3E = -A = -I3F (178) 

and S = {(3/2)M where M is the generalization of the Komar mass, defined as the 
thermally weighted integral 

f3M = 2 [ (3ioceVhd'^-^x (179) 



All these results have been obtained previously in the context of general relativity in 
ref.[69] and the above analysis shows that they continue to remain valid for Lanczos- 
Lovelock models, strengthening the thermodynamic connection. In the context of 
general relativity, using S = Ah/^L'\, where Ah is the horizon area, the relation 
5" = {(3/2)M can be written as M = {l/2)nkBT with n = An/Lj:,. This is just the 
law of equipartition of energy among the horizon degrees of freedom [131] if we think 
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of horizon as divided into patches of size Lp. More generally, the definition of entropy 
given in ref.[69] leads to the equipartition law as an integral over the local acceleration 
temperature Tioc = {Na^n^) /2tt and is given by 



thereby identifying the number of degrees of freedom to be dn = y/a d'^x/ L"], in an area 
element y/ad^x. (One can, alternatively, write down this relation by inspection and 
obtain the gravitational field equations as a consequence, provided one accepts the choice 
of various numerical factors.) This suggests that we interpret the relation S = {/3/2)M 
as law of equipartition even in the more general context of Lanczos-Lovelock theories 
and identify M / (^kBT) = AS as the effective number of degrees of freedom. (In the 
Lanczos-Lovelock models this will not be proportional to the area.) The Noether current 
and its relation to entropy will play a crucial role in our future discussions. 

5.9. Why are gravitational actions Junctionals holographic? 

We shall now reinterpret the above results by approaching them, from first principles, 
in a manner which makes these relations logically transparent [39] . 

We will start with the principle of equivalence and draw from it three important 
consequences. First, it implies that, at least in the long wavelength limit, gravitational 
field has a geometrical description in terms of the metric of the spacetime [132]. Second, 
the validity of laws of special relativity in the local inertial frames allows one to determine 
the influence of gravity on light rays. From this we conclude that, in any given spacetime, 
there will exist observers who will have access to only part of the spacetime. That 
is, horizons — as perceived by a class of observers — will exist in any geometrical 
description of gravity that obeys principle of equivalence. Finally, the principle of 
equivalence also demands (indirectly) a principle of 'democracy of observers'. In flat 
spacetime, we can choose a special coordinate system with the global metric being rjab] 
so when the spacetime metric is given to be gab{t,^) we can always attribute the part 
{dab — Vab) to the choice of non-inertial coordinates. In a curved spacetime we cannot 
do this and hence it no longer makes sense to ask "how much of gat" is due to our using 
non-inertial coordinates and "how much" is due to genuine gravity. Different observers 
in different states of motion might use different coordinates leading to different sets of 
ten functions for gab(t,x.). Because we have no absolute reference metric, it follows that 
no coordinate system or observer is special and the laws of physics should not select out 
any special class of observers. 

Therefore, we need a principle to handle the fact that different observers will have 
access to different regions of the spacetime [133]. If a class of observers perceive a 
horizon, they should still be able to do physics using only the variables accessible to 
them without having to know what happens on the other side of the horizon. Classically, 
the light cone structure dictated by the horizon ensures that the region inside the horizon 
cannot influence the outside. This is, however, not true quantum mechanically when, 




(180) 
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for example, we take into account the entanglement of fields across the horizon. Hence 
there should exist a mechanism which will encode the information in the region V, which 
is inaccessible to a particular observer, at the boundary dV of that region. 

One possible way of ensuring this is to have a suitable boundary term in the action 
principle which will contain the necessary information for observers who perceive a 
horizon. Such a possibility leads to three immediate consequences. First, if the theory 
is generally covariant, so that observers with horizons (like, for example, uniformly 
accelerated observers using a Rindler metric) need to be accommodated in the theory, 
such a theory must have an action functional that contains a surface term. (Since 
general covariance leads to the Noether current, which in turn shows that the action 
will have surface term — see Eq. (172) — we already have a direct demonstration of 
this connection in all static spacetimes.) Second, if the surface term has to encode the 
information which is blocked by the horizon, then there must exist a simple relation 
between the bulk term and surface term in the action and it cannot be arbitrary. Third, 
if the surface term encodes information which is blocked by the horizon, then it should 
actually lead to the entropy of the horizon. In other words, we should be able to compute 
the horizon entropy by evaluating the surface term. 

These three requirements are very strong constraints on the nature of action 
functional describing a theory of gravity. In fact these are sufficiently powerful for 
us to reconstruct the action functional for this class of theories. We shall illustrate 
[134, 130] how this can be done in the case of Einstein's theory. 

To do this we shall start with the assumption that the action in general 
relativity will have a surface term (obtained by integrating a local divergence) which is 
holographically related to the bulk term by Eq. (148). This specific form is needed to 
ensure that the same equations of motion are obtained from Abuik ot from another A' 
(both, as yet, unknown) where: 



A' = d a;A/^Lbuik - d xd, 



9ab- 



d{dcgab) 

Abulk + / d^xd^iV^V'') (181) 



with constructed from gab and F*;, with no further explicit dependence on y/^, 
which has been factored out. Further, must be linear in the F's since the original 
Lagrangian Lbuik was quadratic in the first derivatives of the metric. Since Fs vanish 
in the local inertial frame and the metric reduces to the Lorentzian form, the action 
^buik cannot be generally covariant. Our aim is to determine the surface term using 
the known expression for horizon entropy and then determine the bulk term which is 
consistent with the holographic relation. At this stage we have no assurance that the 
resulting action A' will be generally covariant but it is an important consistency check 
on the idea. 

To obtain a quantity V^, which is linear in Fs and having a single index c, from gab 
and T^ji^, we must contract on two of the indices on F using the metric tensor. (Note 
that we require Abuik, A' etc. to be Lorentz scalars and P'^, V"^ etc. to be vectors under 
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Lorentz transformation.) Hence the most general choice for is the hnear combination 



[(^l9 J- km 



+ Ci29''^k) (182) 

where ai and 02 are two constants. Using the identities = c?fc(ln y/^), y/—gg^^T1f^ = 
—db{y/—gg'"^), we can rewrite the expression for P'^ = ^—gV^ as 

p- = ^V"^ = c,g'%V^ + C2V^dtg''' (183) 

where Ci = Oi — 02, C2 = — 02 are two other unknown constants. We now fix these 
coefficients by using our demand that A^^,. should lead to the entropy of the horizon, 
thereby determining the surface term and — by integrating — the Lagrangian Lbuik- 
We will use a Euclidean Rindler frame for this purpose. 

The form of Rindler metric that is sufficiently general for our purpose can be taken 
to be 



ds^ = 2kI dt%-\ + (dy^ + dz"^) 



2kI 



:i84) 



2kI{x) dt\ + 



Pdx' 



{dy"^ + dz^ 



2kI{x) 

where l{x) is an arbitrary function and /' = (dl/dx). The metric in the first line is same 
as the one in Eq. (7); the second line introduces an arbitrary coordinate x through the 
function l{x). The horizon is at l{x) = and k is the surface gravity for all l{x), giving 
the horizon temperature to be K,/27r. Evaluating the surface term P^ in Eq. (183) for 
this metric, we get the only nonzero component to be 



2k 



C2 



2co 



so that the surface term in the action in Eq. (181) becomes 



A. 



/3P^ / d^xi_ = PP^Aa 



4:7lAl 



C2 



;i85) 



:i86) 



where A± is the transverse area of the {y — z) plane and the time integration is limited to 
the range (0, /3) with /3 = 211/1%. Demanding that A^^^ should give the horizon entropy, 
which we take to be proportional to the transverse area, as —{A^/4:Ap) where Ap is a 
constant, we get the condition: 

1 



C2 + -p-[Ci 



2C2 



;i87) 



IQtxAp 

(The minus sign arises because we are working with Euclidean signature as in, for 
example, Eq. (152).) We demand that this equation should hold, on the horizon, for 
all functions l{x). If we take l{x) = x", then II" /V^ = {n — l)/n. Taking n = 1 makes 
the second term on the left hand side of Eq. (187) vanish giving C2 = — (167r^p)~^. For 
other values of n, we now require the second term to vanish identically which is possible 
only if Ci = 2c2. (One might have thought that, since the horizon is at / = 0, the second 
term on the left hand side of Eq. (187) vanishes identically on the horizon and we can 
only determine C2. The above analysis shows that, for l{x) = with n > 0, the horizon 
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is still at / = X = but II" /V^ = {n — \)/n does not vanish on the horizon. It is this fact 
which allows us to determine Ci and C2. In fact we only need to use n = 1 and n = 2 
corresponding to the two standard forms of the Rindler metric to fix the two constants.) 
This completely determines Ci and C2 and hence the surface term. Continuing back to 
Lorentzian sector, we get 

^ ^ -Mgg'l (188) 



IQnAp \f--g 

This is precisely the surface term in Einstein-Hilbert action as can be seen by comparing 
with Eq. (147) and recalling P'^ = y/—g V^. Given the form of P'^ we need to solve the 
equation 



9ab = = , ^ W'') (189) 



dgab,c J IGirAp y/^ 

to obtain the first order Lagrangian density. It is straightforward to show that this 
equation is satisfied by the Lagrangian 

= (v^ g'' {TTiTL - rf.r- )) . (i90) 

which we already know from Eq. (148). (The solution to Eq. (189) obtained in Eq. (190) 
is not unique. However, self consistency requires that the final equations of motion for 
gravity must admit the line element in Eq. (184) as a solution. It can be shown, by 
fairly detailed algebra, that this condition makes the Lagrangian in Eq. (190) to be 
the only solution.) Given the two pieces, the final second order Lagrangian follows 
from our Eq. (181) and is, of course, the standard Einstein-Hilbert Lagrangian. In this 
approach, our full second order Lagrangian turns out to be the standard Einstein-Hilbert 
Lagrangian. Our result has been obtained, by relating the surface term in the action to 
the entropy per unit area; i.e., the surface terms dictates the form of the Lagrangian in 
the bulk through the holographic relation. 

The same procedure works for Lanczos-Lovelock models though the mathematics 
is much more tedious. (This approach has also uncovered several other issues related to 
entropy, quasi-normal modes etc. and even a possibility of entropy being quantized [135] 
but we will not discuss these aspects.) Also, in the case of Lanczos-Lovelock gravity, 
the expression for entropy has no simple physical motivation unlike in Einstein's theory. 
Hence, while this does illustrate the power of the holographic action principle it does 
not allow one to make significant further progress. We shall see later (see Sec 7.3) that 
the thermodynamic interpretation of the field equations actually arises from different 
approach in which we do not consider the metric as a fundamental variable. 



5.10. Summary 



The various results obtained in the previous sections can be briefly summarized as 
follows. 
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(a) It is possible to associate the notion of temperature with any bifurcation horizon 
in a fairly general manner by, for example, using the periodicity in the imaginary time. 
Association of entropy with a generic horizon is conceptually more involved but it seems 
quite unnatural to associate temperature with all horizons and entropy with only a 
subset of them. The association of thermodynamic variables to a horizon is not a 
special feature of Einstein's theory and extends to a much wider class of theories of 
gravity. 

(b) There are strong hints which suggests that horizon thermodynamics has a 
deep connection with gravitational dynamics, which is not apparent in the conventional 
approach to gravity. In particular, we have seen that: 

• The field equations of gravity in a very wide class of theories reduce to the 
thermodynamic identity TdS = dE + PdV on the horizon. This result has no 
explanation in the conventional approach. 

• The action functional in a wide class of theories of gravity contains both a bulk 
term and a surface term. There is a specific relationship between these two terms 
which allows these action functionals to be interpreted as a momentum space action 
functionals. It is not clear why this peculiar relation exists between the bulk and 
surface terms. 

• The surface term of the action functional, when evaluated on the horizon in a 
solution, gives the entropy. This is quite mysterious since the field equations are 
obtained by varying the bulk term after the surface term has been ignored (or 
cancelled out by a counter-term). Therefore, the field equations and their solutions 
are completely independent of the surface term. We do not expect a specific 
property of the solution (for example, the horizon entropy) to be obtainable from 
the surface term which played no role in the field equations. 

• The Euclidean action in any static spacetime can be interpreted as the free energy 
in a wide class of theories of gravity showing that the minimization of the action 
can be related to the minimization of free energy. 

In the next two sections of this review, we shall provide an alternate perspective 
on gravity which will help us to understand these features better and in a fairly unified 
manner. 

6. The emergent spacetime 

In this section and the next, we shall present an alternative perspective on the nature 
of gravity motivated by the results described in the previous section. This approach, as 
we shall see, provides a natural setting for many of the results obtained in the earlier 
sections which are somewhat mysterious in the context of the conventional approach. 
The key feature of the new perspective is that it treats gravity as an emergent phenomena 
and describes its dynamics in the thermodynamic limit. We will, therefore, begin by 
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making clear what is meant by emergent phenomena in this context and estabhshing its 
connection with thermodynamics. 

As we have already argued in Sec. 1, the fact that spacetimes can be hot strongly 
indicates the existence of internal degrees of freedom for the spacetime. This fact 
probably would be accepted by most people working in quantum gravity since almost all 
these models introduces extra structures at microscopic (Planck) scales in the spacetime. 
The natural picture which then emerges is that there are some "atoms of spacetime" at 
the microscopic level and the description of spacetime in terms of variables like metric, 
curvature etc. is a continuum, long wavelength, approximation. This is analogous to 
description of a gas or a fluid in terms of dynamical variables like density p, velocity 
V etc. in the continuum limit, none of which have any relevance in the microscopic 
description. 

The new ingredient we will introduce is based on the fact that, while we do not have 
definite knowledge about the statistical mechanics of atoms of spacetime, we should 
be able to develop the thermodynamic limit of the theory taking clues from horizon 
thermodynamics. 

As emphasized in Sec. 1, such an approach has two major advantages. First, 
the thermodynamic description has a universal validity which is fairly independent of 
the actual nature of the microscopic degrees of freedom. Second, the entropy of the 
system arises due to our ignoring the microscopic degrees of freedom. Turning this 
around, one can expect the form of entropy functional to encode the essential aspects 
of microscopic degrees of freedom, even if we do not know what they are. If we can 
arrive at the appropriate form of entropy functional, in terms of some effective degrees 
of freedom using our knowledge of horizon thermodynamics, then we can expect it to 
provide the correct description.^ As we know, thermodynamics was developed and 
used effectively decades before we understood the molecular structure of matter or its 
statistical mechanics. 

Similarly, even without knowing the microstructure of spacetime or the full quantum 
theory of gravity, we should be able to make significant progress with the thermodynamic 
description of spacetime. The horizon thermodynamics provides [70] valuable insights 
about the nature of gravity totally independent of what "the atoms of spacetime" may 
be. This is what we will attempt to do. (There has been several other attempts 
in literature to implement the idea that gravity is an emergent phenomenon, which 
we shall not discuss. They do not: (a) address issues we have raised regarding the 
action functionals and (b) cannot handle Lanczos-Lovelock models effectively; for a 
small sample of papers on other approaches, which contain additional references, see 
Ref.[136].) 

% Incidentally, this is why thermodynamics needed no modification due to either relativity or quantum 
theory. An equation like TdS = dE+PdV will have universal applicability as long as effects of relativity 
or quantum theory are incorporated in the definition of S{E, V) appropriately. 
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6.1. Thermodynamic interpretation of field equations of gravity 

Consider the action functional in Eq. (79) which, on variation, leads to the field equations 

26^ab-T,, = (191) 

where the explicit form of Qab is given by Eq. (81). As mentioned earlier, this equation 
— which equates a geometrical quantity to matter variables — does not have any simple 
physical interpretation. The lack of an elegant principle to determine the dynamics of 
gravity is in sharp contrast with the issue of determining the kinematics of gravity which 
can be tackled through the principle of equivalence by demanding that all freely falling 
observers must find that the equations of motion for matter degrees of matter must 
reduce to their special relativistic form. 

Our first aim will be to remedy this situation and provide a physical interpretation 
to Eq. (191). We will do this in a manner very similar in spirit to using freely 
falling observers to determine the kinematics of gravity. At every event in spacetime, 
we will introduce uniformly accelerating local Rindler observers and use the horizon 
thermodynamics perceived by these Rindler observers to constrain the background 
geometry. We shall begin by making the notion of local Rindler observers and their 
coordinate system well defined. 

Let us choose any event V and introduce a local inertial frame (LIE) around it with 
Riemann normal coordinates X" = (T, X) such that V has the coordinates X"- = in 
the LIE. Let k"" be a future directed null vector at V and we align the coordinates of 
LIE such that it lies in the X — T plane at V. We next transform from the LIE to 
a local Rindler frame (LRE) coordinates by accelerating along the X-axis with an 
acceleration k by the usual transformation. The metric near the origin now reduces to 
the form 

ds"^ = -dT'^ + dX^ + dxl = -K^x^dt^ + dx"^ + dLy.\ 
df 

= - 2kI dt^ + — + c/xi (192) 
2 m 

where (t, /,xx) and (t, x,xx) are the coordinates of LRE. Let be the approximate 
Killing vector corresponding to translation in the Rindler time such that the vanishing 
of C"^a = —N"^ characterizes the location of the local horizon "H in LRE. Usually, we 
shall do all the computation on a timelike surface infinitesimally away from with = 
constant, called a "stretched horizon". Let the timelike unit normal to the stretched 
horizon be r^. 

This LRE (with metric in Eq. (192)) and its local horizon "H will exist within 
a region of size L < 7^~l/2 (where 7^ is a typical component of curvature tensor of 
the background spacetime) as long as <^ 'Rr^l'^ . This condition can always be 
satisfied by taking a sufficiently large k (see Eig. 3). This procedure introduces a class 
of uniformly accelerated observers who will perceive the null surface T = ±X as the 
local Rindler horizon "H. 

Essentially, the introduction of the LRE uses the fact that we have two length 
scales in the problem at any event. Eirst is the length scale TZ~^^'^ associated with the 
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Figure 3. The left frame illustrates schematically the light rays near an event V in the 
t — X plane of an arbitrary spacetime. The right frame shows the same neighbourhood 
of V in the locally inertial frame at V in Riemann normal coordinates (T^X). The 
light rays now become 45 degree lines and the trajectory of the local Rindler observer 
becomes a hyperbola very close to T = zLX lines which act as a local horizon to the 
Rindler observer. 



Figure 4. The region around V shown in Fig. 3 is represented in the Euclidean sector 
obtained by analytically continuing to imaginary values of T by Te = iT. The horizons 
T = ±X collapse to the origin and the hyperbolic trajectory of the Rindler observer 
becomes a circle of radius around the origin. The Rindler coordinates {t, x) become 
— on analytic continuation to t^; = it — the polar coordinates (r = x,6 — kIe) near 
the origin. 



curvature components of the background metric over which we have no control while 
the second is the length scale associated with the accelerated trajectory which we 
can choose. Hence we can always ensure that <^ IZ^^^"^. In fact, this is clearly 
seen in the Euclidean sector in which the horizon maps to the origin (see Fig. 4). The 
locally flat frame in the Euclidean sector will exist in a region of radius TZ^^^'^ while the 
trajectory of a uniformly accelerated observer will be a circle of radius k,~^ and hence 
one can always keep the latter inside the former. The metric in Eq. (192) is just the 
metric of the locally flat region in polar coordinates. 

More generally, one can choose a trajectory x*(r) such that its acceleration 
= u'^ViU^ (where is the time-like four velocity) satisfies the condition a^aj = k^. 
In a suitably chosen LIE this trajectory will reduce to the standard hyperbola of a 
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uniformly accelerated observer. 

Our construction also defines local Rindler horizons around any event. Further, 
the local temperature on the stretched horizon will be k/2tiN so that /3ioc = PN with 
/3 = h/2tt. Note that in the Euclidean sector the Rindler observer's trajectory is a circle 
of radius which can be made arbitrarily close to the origin. Suppose the observer's 
trajectory has the usual form X = k,~^ cosh nt; T = sinh nt which is maintained for 
a time interval of the order of 271 /n. Then, the trajectory will complete a full circle in 
the Euclidean sector irrespective of what happens later. When we work in the limit of 
K — oo, our construction becomes arbitrarily local in both space and time [60, 137]. 

As stressed earlier in Sec. 4.4, the local Rindler observers will perceive the 
thermodynamics of matter around them very differently from the freely falling observers. 
In particular, they will attribute a loss of entropy SS = {2'k/h)5E (see Eq. (113)) 
when an amount of energy 5E gets close to the horizon (within a few Planck lengths, 
say). In the Rindler frame the appropriate energy- momentum density is T^^'^- (It 
is the integral of T^^'^d'La that gives the Rindler Hamiltonian Hr, which leads to 
evolution in Rindler time t and appears in the thermal density matrix p = exp —(5Hr.) 
A local Rindler observer, moving along the orbits of the Killing vector field with 
four velocity = ^""/N, will associate an energy density u"'{TabC,^) and an energy 
6E = u'^(Tab^^)dVpj-op with a proper volume dVpj-op- If this energy gets transfered 
across the horizon, the corresponding entropy transfer will be 55'matter = PiocSE where 
Aoc = /3A^ = {27t/k)N is the local (redshifted) temperature of the horizon and is the 
lapse function. Since AocM" = il3N){^^/N) = (3^", we find that 

55mattcr = (^C^^Tab dVp.op (193) 

Consider now the gravitational entropy associated with the local horizon. From the 
discussion of Noether charge as horizon entropy in Section 3.5 [see Eq. (93)], we know 
that f3iocJ°', associated with the Killing vector can be thought of as local entropy 
current. Therefore, SS = f3iocUaJ°'dVprop can be interpreted as the gravitational entropy 
associated with a volume dVprop as measured by an observer with four-velocity m". 
(The conservation of J" ensures that there is no irreversible entropy production in 
the spacetime.) Since is a Killing vector locally, satisfying Eq. (88) it follows that 
6^v = giving the current to be J" = [L^"" + 2Q°-^^h}- For observers moving along the 
orbits of the Killing vector with u"' = C,"'/N we get 

55grav = PiocNUardVp^^p = PWdVp^.p = l5[i^U'^g'^')+L{i,i')] f/V^prop(194) 

As one approaches the horizon, ^"^^ — >■ making the second term vanish and we find 
that 

5^grav = /3[er(26^a,)] dVp.^p (195) 

In the same limit will become proportional to the original null vector y we started 
with keeping everything finite. We now see that the condition 5S^ra.w = SSuiatter leads to 
the result 

[2gab-TabWk' = (196) 
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Since the original null vector ka was arbitrary, this equation should hold for all null 
vectors for all events in the spacetime. This is equivalent to 2Q°-'^ — T°'^ = Xg"'^ with 
some constant A. (Because of the conditions Va^"'' = 0, VaT"^^ = 0, it follows that A 
must be a constant.) 

We have thus succeeded in providing a purely thermodynamical interpretation 
of the field equations of any diffeomorphism invariant theory of gravity. Note that 
the equations Eq. (196) has an extra symmetry which standard gravitational field 
equations do not have: This equation is invariant under the shift T"^ — > T"^ + ^g"-^ 
with some constant This symmetry has important implications for cosmological 
constant problem which we will discuss in Section 7.5. 

It is clear that the properties of LRF are relevant conceptually to define the 
intermediate notions (local Killing vector, horizon temperature ....) but the essential 
result is independent of these notions. Just as we introduce local inertial frame to decide 
how gravity couples to matter, we use local Rindler frames to interpret the physical 
content of the field equations. 

In Section 3.5 we mentioned that J" is not unique and one can add to it the 
divergence of any anti-symmetric tensor (see the discussion after Eq. (85)). In providing 
the thermodynamic interpretation to the field equations, we have ignored this ambiguity 
and used the expression in Eq. (83). There are several reasons why the ambiguity is 
irrelevant for our purpose. First, in a truly thermodynamic approach, one specifies 
the system by specifying a thermodynamic potential, say, the entropy functional. In a 
local description, this translates into specifying the entropy current which determines 
the theory. So it is perfectly acceptable to make a specific choice for consistent 
with the symmetries of the problem. Second, we shall often be interested in theories 
in which the equations of motion are no higher than second order and has the form 
in Eq. (97). In these Lanczos-Lovelock models it is not natural to add any extra term 
to the Noether current such that it is linear in as we approach the horizon with a 
coefficient determined entirely from metric and curvature. Finally, we shall obtain in 
Section 7 the field equation from maximizing an entropy functional when this ambiguity 
will not arise. 

It may be noted that our result only required the part of J" given by 2^fe"^^ Hence, 
one can obtain the same results by postulating the entropy current to be J*^ = 2Q^^^ 
which is also conserved off-shell when is a Killing vector. In fact one can give 2Q^^^ 
an interesting interpretation. Suppose there are some microscopic degrees of freedom 
in spacetime, just as there are atoms in a solid. If one considers an elastic deformation 
— x° + ^"(x) of the solid, the physics can be formulated in terms of the displacement 
field ^"(x) and one can ask how thermodynamic potentials like entropy change under 
such displacement. Similarly, in the case of spacetime, one could think of 

^Sgrav = (3loc{'2Gt)Ua6x'' (197) 

as the change in the gravitational entropy density under the 'deformation' of the 
spacetime a;" — + dx"" as measured by the Rindler observer with velocity u". (One 



Thermodynamical Aspects of Gravity: New insights 



65 



can show that this interpretation is consistent with all that we know about horizon 
thermodynamics.) So, one can interpret the left hand side of gravitational field equation 
(2Q^) as giving the response of the spacetime entropy to the deformations. This matches 
with the previous results because Aoc^a = /3^a implies that the entropy density will be 
proportional to 2Qabk°'k^ on the horizon. We will see later that this interpretation 
remains valid in a very general context. 

Once it is understood that the real physical meaning of the field equations lies in 
writing them in the form of Eq. (196), it is possible to re-interpret these equations in 
several alternative ways all of which have the same physical content. We shall mention 
two of them. 

Consider an observer who sees some matter energy flux crossing the horizon. Let 
Ta be the spacelike unit normal to the stretched horizon E, pointing in the direction of 
increasing A^. The energy flux through a patch of stretched horizon will be Tab^°'r^ 
and the associated entropy flux will be PiocTabi""^^ ■ To maintain the second law 
of thermodynamics, this entropy flux must match the entropy change of the locally 
perceived horizon. The gravitational entropy current is given by PiocJ"'-, such that 
PiociTaJ"") gives the corresponding gravitational entropy flux. So we require 



to hold at all events where is the conserved Noether current corresponding to 
The product r^J" for the vector r", which satisfies ^""Ta = on the stretched horizon is 
r^J" = 2Q"-^raib- Hence we get 



As -^■ and the stretched horizon approaches the local horizon and A^r* approaches 
^* (which in turn is proportional to /c*) so that Pioc^a = PNva Pia- So, as we approach 
the horizon Eq. (199) reduces to Eq. (196). 

There is another way of interpreting this result which will be useful for further 
generalizations. Instead of allowing matter to flow across the horizon, one can equally 
well consider a virtual, infinitesimal (Planck scale), displacement of the "H normal to 
itself engulfing some matter. We only need to consider infinitesimal displacements 
because the entropy of the matter is not 'lost' until it crosses the horizon; that is, 
until when the matter is at an infinitesimal distance (a few Planck lengths) from the 
horizon. Some entropy will be again lost to the outside observers unless displacing a 
piece of local Rindler horizon costs some entropy. 

We can verify this as follows: An infinitesimal displacement of a local patch of the 
stretched horizon in the direction of r^, by an infinitesimal proper distance e, will change 
the proper volume by dVprop = ey/ad^~'^x where aab is the metric in the transverse space. 
The flux of energy through the surface will be T^^^^Va and the corresponding entropy 
flux can be obtained by multiplying the energy flux by /3ioc. Hence the 'loss' of matter 
entropy to the outside observer because the virtual displacement of the horizon has 
engulfed some matter is 5Sm = Pxoc^E = PiocT""^ jdVprop- To find the change in the 



(198) 



(199) 
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gravitational entropy, we again use the Noether current J" corresponding to the local 
Killing vector Multiplying by r" and /3ioc = PN, we get 

Accra J" = PiocUaT''' + l3N{raC)L (200) 

As the stretched horizon approaches the true horizon, we know that Nr"" — )■ and 
l3^°'^aL — )■ making the last term vanish. So 

'^'S'grav = f^^aJ dVprop = f3T ^a^jdVprop = (201) 

showing the validity of local entropy balance for any (5. In this limit, ^* also goes to 
kA/c* where A is the affine parameter associated with the null vector k"' we started with 
and all the reference to LRF goes away. 

7. Gravity: The inside story 

7.1. An Entropy maximization principle for gravitational field equations 

The last interpretation given above is similar to switching from a passive point of view to 
an active point of view. Instead of allowing matter to fall into the horizon, we are making 
a displacement of the horizon surface to engulf the matter when it is infinitesimally close 
to the horizon. But in the process, we have introduced the notion of virtual displacement 
of horizons and, for the theory to be consistent, this displacement of these surface degrees 
of freedom should cost some entropy. This allows one to associate an entropy functional 
with the normal displacement of any horizon. 

An analogy may be helpful in this context. If gravity is an emergent, long 
wavelength, phenomenon like elasticity then the diffeomorphism x°' — + ^" is 
analogous to the elastic deformations of the "spacetime solid" [138]. It then makes sense 
to demand that the entropy density should be a functional of and its derivatives V^^". 
By constraining the functional form of this entropy density, we should be able to obtain 
the field equations of gravity by a maximization principle. Recall that thermodynamics 
relies entirely on the form of the entropy functional to make predictions. Hence, if we 
can determine the form of entropy functional for gravity (Sgrav) in terms of the normal to 
the null surface, then it seems natural to demand that the dynamics should follow from 
the extremum prescription S[Sgrav + Smatter] = for all null surfaces in the spacetime 
where Smatter is the relevant matter entropy. 

The form of Smatter and S'g^av can be determined as follows. Let us begin with 
Smatter which is casy to ascertain from the previous discussion. If Tab is the matter 
energy- momentum tensor in a general D{> A) dimensional spacetime then an expression 
for matter entropy relevant for our purpose can be taken to be 

5matt = / d^'x^Tabu'^n' (202) 
Jv 

where is a null vector field. From our Eq. (193) we see that the entropy density 
associated with proper 3- volume is f3{Tab^°'^'')dVprop where — on the horizon — the 
vector becomes proportional to a null vector n". If we now use the Rindler coordinates 



Thermodynamical Aspects of Gravity: New insights 



67 



in Eq. (7) in which y/—g = 1 and interpret the factor /3 as arising from an integration 
of dt in the range (0,/9) we find that the entropy density associated with a proper 
four volume is {Tahn°"n^)- This suggests treating Eq. (202) as the matter entropy. For 
example, if Tab is due to an ideal fluid at rest in the LIE then Tabn°"n^ will contribute 
{p + P), which — by Gibbs-Duhem relation — is just Tiocais where s is the entropy 
density and TJ~2ii ~ /^-^ properly redshifted temperature. Then 



which matches with Eq. (202) in the appropriate limit. 

It should be stressed that this argument works for any matter source, not necessarily 
the ones with which we conventionally associate an entropy. What is really relevant is 
only the energy flux close to the horizon from which one can obtain an entropy flux. 
We do have the notion of energy flux across a surface with normal being Tah^^r"" 
which holds for any source T"*. Given some energy flux 5E in the Rindler frame, there 
is an associated entropy flux loss 5S = (55E as given by Eq. (113). (One might think, 
at first sight, that an ordered field, say, a scalar field, has no temperature or entropy 
but a Rindler observer will say something different. For any state, she will have a 
corresponding density matrix p and an entropy — Tr(p In p); after all, she will attribute 
entropy even to vacuum state.) It is this entropy which is given by Eq. (113) and 
Eq. (202). The only non-trivial feature in Eq. (202) is the integration range for time 
which is limited to (0,/3). This is done by considering the integrals in the Euclidean 
sector and rotating back to the Lorentzian sector but the same result can be obtained 
working entirely in the Euclidean sector. (There is an ambiguity in the overall scaling 
of since if is null so is f{x)n°' for all /(x); we will comment on this ambiguity, 
which anyway turns out to be irrelevant, later on.) 

Next, let us consider the expression for S'grav We will first describe the simplest 
possible choice and then consider a more general expression. The simplest choice is to 
postulate Sgrav to be a quadratic expression [139] in the derivatives of the normal: 



Jv 

where the explicit form of P^h'^'^ is ascertained below. The expression for the total 
entropy now becomes: 



We should be able to determine the field equations of gravity by extremizing this 
entropy functional. However, there is one crucial conceptual difference between the 
extremum principle introduced here and the conventional one. Usually, given a set of 
dynamical variables Ua and a functional S'[na], the extremum principle will give a set 
of equations for the dynamical variable Ua- Here the situation is completely different. 





(204) 




(205) 
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We expect the variational principle to hold for all null vectors n"" thereby leading to 
a condition on the background metric. Obviously, the functional in Eq. (205) must be 
rather special to accomplish this and one needs to impose restrictions on P^i,'^'^ — and 
Tab though that condition turns out to be trivial — to achieve this. (Of course, one can 
specify any null vector ?t,"(x) by giving its components f^{x) = n'^e^ with respect to 
fixed set of basis vectors with e^ef = 5^ etc so that = f^e\. So the class of all 
null vectors can be mapped to the scalar functions with the condition fAf^ = 0.) 

It turns out — as we shall see below — that two conditions are sufficient to ensure 
this. First, the tensor Pabcd should have the same algebraic symmetries as the Riemann 
tensor Rabcd of the D- dimensional spacetime. This condition can be ensured if we define 
Pj"'^ as 

8T 

Pa'"' = ^ (206) 

where L = L(i?'^^^^, g'^'') is some scalar. The motivation for this choice arises from the fact 
that this approach leads to the same field equations as the one with L as gravitational 
Lagrangian in the conventional approach (which explains the choice of the symbol L). 
Second, we will postulate the condition: 

^^pabcd ^ Q ^207) 

as well as VqT'^^ = which is anyway satisfied by any matter energy-momentum tensor. 

One possible motivation for this condition in Eq. (207) arises from the following fact: 
It will ensure that the field equations do not contain any derivative of the metric which 
is of higher order than second. Another possible interpretation arises from the analogy 
introduced earlier. If we think of n"" as analogous to deformation field in elasticity, 
then, in theory of elasticity [140] one usually postulates the form of the thermodynamic 
potentials which are quadratic in first derivatives of Ua- The coefficients of this term will 
be the elastic constants. Here the coefficients are P"*'^'^ and the condition in Eq. (207) 
may be interpreted as saying the 'elastic constants of spacetime solid' are actually 
'constants'. This is, however, not a crucial condition and in fact we will see below 
how this condition in Eq. (207) can be relaxed. 



7.2. The field equations 

Varying the normal vector field n°- in Eq. (205) after adding a Lagrange multiplier 
function A(x) for imposing the condition UaSn"' = 0, we get 

-SS = 2 [ d^x^ [APab^'^V ji" {Vd5n^) - Tab^fdn!' 
Jv 

-\{x)gabn''5n^] (208) 

where we have used the symmetries of P^b'^'^ and Tab- (We note, for future reference, 
that the Lagrange multiplier in the calculation only imposes the constancy of n^n* under 
variation and does not require rij to be null vector.) An integration by parts and the 
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condition Vd.Pab^'^ — 0' leads to 

-6S = 2 [ rf^xv^ [-^Pab"' (VdV.n'^) - (T,, + XgabK] Sn' 
Jv 

+ 8 [ d'^~^xVh[kdPab'"^{Vcn'')]5n\ (209) 

JdV 

where /c" is the /^-vector field normal to the boundary dV and h is the determinant 
of the induced metric on dV. As usual, in order for the variational principle to be 
well defined, we require that the variation 6n°' of the vector field should vanish on the 
boundary. The second term in Eq. (209) therefore vanishes, and the condition that 
S[n°-] be an extremum for arbitrary variations of then becomes 

2P,,^^ (VeVd - VrfVj - (T,, + Xg,,)n'' = , (210) 

where we used the antisymmetry of P^b'^'^ in its upper two indices to write the first term. 
Using the definition of the Riemann tensor in terms of the commutator of covariant 
derivatives and writing TZ^ = Pb^^R^iju the above expression reduces to 

{2ni -Tl + \5l) = {2QI - T," + (L + = , (211) 

where we have used the definition of Ql in Eq. (97). We see that the equations of motion 
do not contain derivatives with respect to n°- which is, of course, the crucial point. This 
peculiar feature arose because of the symmetry requirements we imposed on the tensor 
Pab'^'^- (Multiplying by and noting n°'na = we see that Eq. (211) and Eq. (196) 
are identical.) We need the condition in Eq. (211) holds for arbitrary vector fields n". 
One can easily show [139] using VaGb = = VaT^" that this requires A + L = constant 
leading to the field equation 



Gb 



K - l^bL 



\Tt + (212) 



where A is a constant. Comparison of Eq. (97) (or Eq. (81)) with Eq. (212) shows 
that these are precisely the field equations for gravity in a theory with Lagrangian L 
when Eq. (205) is satisfied. One crucial difference between the two equations is the 
introduction of the cosmological constant A as an integration constant in Eq. (212); we 
will discuss this later in Section 7.5. 

We mentioned earlier that the expression in Eq. (205) depends on the overall scaling 
of which is arbitrary, since f{x)n°' is a null vector if n'^ is null. But since the arbitrary 
variation of with the constraint Uau"- = includes scaling variations of the type 
(5n" = e(x)n", it is clear that this ambiguity is irrelevant for determining the equations 
of motion. 

To summarize, we have proved the following. Suppose we start with a total 
Lagrangian L{Rabcd, Qab) + Lmatt, define a P"-^^'^ by Eq. (206) ensuring it satisfies 
Eq. (207). Varying the metric with this action will lead to certain field equations. 
We have now shown that we will get the same field equations (but with a cosmological 
constant) if we start with the expression in Eq. (205), maximize it with respect to 
and demand that it holds for all n". 
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This result might appear a httle mysterious at first sight, but the following 
alternative description will make clear why this works. Note that, using the constraints 
on P"-'"^'^ we can prove the identity 

= 4Ve[P,fe^VVdn''] - 2n-P,,'^Vi,Vd]n'' 

= 4V.[P,,^VV,r2''] - 2n'^P,,^''R\,X 

= 4Vc[Pab"''n''Vdn'] + 2n"e?,,n* (213) 

where the first line uses Eq. (207), the second line uses the antisymmetry of P^b'^'^ in c 
and d, the third line uses the standard identity for commutator of covariant derivatives 
and the last line is based on Eq. (81) when nau"" = and Eq. (207) hold. Using this in 
the expression for S in Eq. (205) and integrating the four-divergence term, we can write 
the entropy functional as 

S[rf]= - f (i^-ixA;,v^(4P^b'=VVdn'') 

JdV 

- [ d''xV^[{2gab-TabKn'] (214) 
Jv 

So, when we consider variations ignoring the surface term we are effectively varying 
{2Qab — Tab)n""n^ with respect to Ua and demanding that it holds for all Ua- This is 
the reason why we get {2Qab = Tab) except for a cosmological constant. There is an 
ambiguity of adding a term of the form \{x)gab in the integrand of the second term in 
Eq. (214) leading to the final equation {2Qab = Tab + \{x)gab) but the Bianchi identity 
V aQ"'^ = along with VqT"'' = will make X{x) actually a constant. (We see from 
Eq. (213) that, in the case of Einstein's theory, we have a bulk Lagrangian n°(V[aVfe])n'' 
for a vector field n° — plus for a surface term which does not contribute to variation. 
In flat spacetime, in which covariant derivatives become partial derivatives, the bulk 
lagrangian becomes vacuous; i.e., there is no bulk dynamics in n", in the usual sense. 
Nevertheless, they do play a crucial role.) 

The expression in Eq. (214) also connects up with our previous use of 2Qabn"'n^ as 
gravitational entropy density. The gravitational part of the entropy in Eq. (214) can be 
written as 

SgraAn""] = - I d^xv^4P,,^'^Ven'^V,n^ (215) 
Jv 

= - I rf^-^xfceV^(4P,,'='^n»Vrfn^) - / rf^xv^(2e;„6n"n'') 

JdV Jv 

with one bulk contribution (proportional to 2Qab'n°'n^) and a surface contribution. When 
equations of motion hold, the bulk also get a contribution from matter which cancels it 
out leaving the entropy of a region V to reside in its boundary dV. 

It is now clear how we can find an S for any theory, even if Eq. (206) does not hold. 
This can be achieved by starting from the expression {2Qab — Tab)n""n^ as the entropy 
density, using Eq. (81) for Qab and integrating by parts. In this case, we get for S'grav 
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the expression: 

Jv 

+ {VcVdP''''"')nanf,] (216) 

Varying this with respect to will then lead to the correct equations of motion and — 
incidentally — the same surface term. 

While one could indeed work with the more general expression in Eq. (216), there 
are four reasons to prefer the imposition of the condition in Eq. (206). First, as we 
shall see below, with that condition we can actually determine the form of L; it turns 
out that in D=4, it uniquely selects Einstein's theory, which is probably a nice feature. 
In higher dimensions, it picks out a very geometrical extension of Einstein's theory 
in the form of Lanczos-Lovelock theories. Second, it is difficult to imagine why the 
terms in Eq. (216) should occur with very specific coefficients. In fact, it is not clear 
why we cannot have derivatives of Rabcd in L, if the derivatives of Pabcd can occur in 
the expression for entropy. Third, it is clear from Eq. (81) that when L depends on 
the curvature tensor and the metric, Qab can depend up to the fourth derivative of the 
metric if Eq. (206) is not satisfied. But when we impose Eq. (206) then we are led to field 
equations which have, at most, second derivatives of the metric tensor which is again a 
desirable feature. Finally, if we take the idea of elastic constants being constants, then 
one is led to Eq. (206). None of these rigorously exclude the possibility in Eq. (216) and 
in fact this model has been explored recently [141]. 

Our variational principle extremises the total entropy of matter and gravity when 
Ua is a null vector. It is, however, possible to provide an alternative interpretation of 
our variational principle (along the lines of ref. [69]), which is of interest when we are 
dealing with static spacetimes with a horizon. Such spacetimes are described by the 
line element 

ds^ = -N\^)dt^ + -i^^{-K)dx^'dx'' (217) 

If n* = ^'^ /N denotes the four velocity of static observers with constant, where 
^* is a timelike Killing vector such that = A^^(x) = is the location of the 
horizon "H, then the matter energy is given by the integral of dU = Tah^""n^ ^d^~^x = 
TabU^n^^f^d^-^x. In this case, our variational principle can be thought of as 
exremising just the gravitational entropy in Eq. (204) subject to two constraints: (i) 
5{nin}) = where n* is now the velocity vector of static observers with n^n* = — 1 and 
(ii) the total matter energy U is constant. Implementing the constancy of U under 
variation by a lagrange multiplier /3 and extremising S — PU, we can as usual identify 
/3 with the range of time integration by analytic continuation from the Euclidean sector 
so that PU becomes an integral over Tabn""n^y/—gd^x. Also note that, when is a 
non-null vector, the identity in Eq. (213) becomes 

AP,,"'Vcn''Vdn' = 4Vc[P,fe'VVdn'] + 2na^n''n' (218) 
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which allows us to work with an alternative definition of 5* given by 

oc j rf^xv^ (27^„^r^V) = /3 j d^-^xVhiJan") (219) 

where the second equality arises on replacing the time integration by multiplication by 
/3 and using ^/—g = Ny/h^n^ = ^^/N along with the expression for Noether current in 
Eq. (89). This result reinforces the idea that this expression is gravitational entropy. In 
the context of Einstein's theory (with TZai = Rai) this reduces to the expression used in 
ref. [69]. More details regarding this approach can be found in ref. [131, 69]. 

Having determined the gravitational field equations, we will make a brief comment 
on the matter sector, before proceeding further. In the conventional action principle, one 
will have a functional which depends on the gravitational degrees of freedom through 
the metric and on the matter degrees of freedom through the matter variables and we 
will vary both to get the equations of motion for gravity and matter. In maximizing the 
entropy we have only varied tIq. However, at the classical level, the equations of motion 
for matter are already contained in the condition VqT"^ = which we have imposed. 
One can do quantum field theory in a curved spacetime using these field equations in 
the Heisenberg picture. Only in the context of path integral quantization of the matter 
fields, one needs to exercise some care. In this case, we should vary n°- first and get 
the classical equations for gravity because the expression in Eq. (205) is designed as an 
entropy functional. But once we have obtained the field equations for gravity, we can 
perform the usual variation of matter Lagrangian in a given curved spacetime and get 
the standard equations [139]. 

So far we have not fixed p°-^'^'^ and so we have not fixed the theory. It is, however, 
possible to determine the form of P"-^'^'^ using Eq. (206) which we shall now describe. 



7.3. The origin of Lanczos- Lovelock models 

In a complete theory, the explicit form of P"^'^'^ will be determined by the long wavelength 
limit of the microscopic theory just as the elastic constants can — in principle — be 
determined from the microscopic theory of the lattice. In the absence of such a theory, we 
need to determine P"^'^'^ by general considerations which is possible when P"-^'^'^ satisfies 
Eq. (207). Since this condition is identically satisfied by Lanczos-Lovelock models which 
are known to be unique, our problem can be completely solved by taking the P"-^'^'^ as a 
series in the powers of derivatives of the metric as: 

P'^'^'^ig^,, R^jki) = ci P "'^"((?.,) + C2 P '^'"="((7.„ R^.ki) + ■ ■ ■ , (220) 

where Ci,C2,-- - are coupling constants with the m th order term derived from the 
Lanczos-Lovelock Lagrangian: 

(ttl) F) T 

^ ah ^ '^abbs...b2m 0314 a2m-ia2»n dR""^ ' 

where ^atbi b2^ alternating tensor. The lowest order term depends only on the 

metric with no derivatives. The next term depends (in addition to metric) linearly on 
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curvature tensor and the next one will be quadratic in curvature etc. The lowest order 
term in Eq. (220) (which leads to Einstein's theory) is 

= ^l^^' = 3^^^^"^^ " ^'^''^ ■ ^^^^^ 

so that when we use Eq. (222) for P^'"''^, Eq. (212) reduces to Einstein's equations. The 
corresponding gravitational entropy functional is: 

SgrIu"] = {Van'Vtn'' - (V^n^f) (223) 

That is, we can obtain the field equations in general relativity by varying the vector 
fields in the above functional and demanding that the resulting equations hold for 
all null vector fields. Interestingly, the integrand in Sgr has the Tr^K"^) — [TrK)'^ 
structure. If we think of the D = 4 spacetime being embedded in a sufficiently large 
k-dimensional flat spacetime we can obtain the same structure using the Gauss-Codazzi 
equations relating the (zero) curvature of k-dimensional space with the curvature of 
spacetime. As mentioned earlier, one can express any vector field in terms of a set 
of basis vector fields e^. Therefore, one can equivalently think of the functional ^gr as 
given by 




(224) 



where P^"^ is a suitable projection operator. It is not clear whether the embedding 
approach leads to any better understanding of the formalism; in particular, it does not 
seem to generalize in a natural fashion to Lanczos-Lovelock models. 

The next order term (which arises from the Gauss-Bonnet Lagrangian) in Eq. (98) 

is: 

= ^ -^^afcaafld T>bzbA 

= ^ [Rf, - G% + G'c^d + Rd^l - (225) 

and similarly for all the higher orders terms. None of them can contribute in Z) = 4 so 
we get Einstein's theory as the unique choice if we assume D = A. If we assume that 
pahcd built only from the metric, then this choice is unique in all D. 



7.4- On-shell value of entropy functional 

The analysis so far used a variational principle based on the functional in Eq. (205). 
While the matter term in this functional has a natural interpretation in terms of entropy 
transfered to the horizon, the interpretation of the gravitational part needs to be made 
explicit. The interpretation of Sgrav as entropy arises from the following two facts. 
First, we see from the identity Eq. (213) that this term differs from 2Qijn^n^ by a total 
divergence. On the other hand, we have seen earlier that the term 2Qijn^n^ can be 
related to the gravitational entropy of the horizon through the Noether current. In fact. 
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Eq. (214) shows that when the equations of motion holds the total entropy of a bulk 
region is entirely on its boundary. Further if we evaluate this boundary term 

-5|on-sheii = 4 / d^'-'xk^Vh {P'^'^^n^VbUd) (226) 

JdV 

(where we have manipulated a few indices using the symmetries of P"*^'^) in the case of 
a stationary horizon which can be locally approximated as Rindler spacetime, one gets 
exactly the Wald entropy of the horizon [139]. 

To prove this, we will use a limiting procedure and provide a physically motivated 
choice of based on the local Rindler frame. Making such a choice is necessary for 
two reasons. First, we do not expect the value of on-shell S to have any direct physical 
meaning for a solution which does not have a horizon. So some choices have to be made. 
Second, we had already mentioned that the expression for S is not invariant under the 
scaling — )■ f{x)n°'. While this is irrelevant for obtaining the field equations, it does 
change the value of on-shell 5*. So we also need to have a prescription for normalization. 
We expect, however, to find sensible results when we evaluate this expression on a local 
Rindler approximation to the horizon which is what we shall do. 

As usual, we shall introduce the LIF and LRF around an event and take the normal 
to the stretched horizon (at = e) to be ra- In the coordinates used in Eq. (192), we 
have the components: 

n„ = (0,1,0,0,...) ; n"^ = (0,1,0,0,...) ; Vh = Ke^/^ , (227) 

where a is the metric determinant of the transverse surface. This vector field n°- is 
a natural choice for evaluation of Eq. (226) if we evaluate the integral on a surface 
with N = e = constant, and take the limit e — )■ at the end of the calculation. In 
the integrand of Eq. (226) for the entropy functional, we use d^~^x = dtd^~^x±, and 
VfeTifi = — F^^^rZa = — F^^^, of which only F^qq = K^e is nonzero. The integrand for the 
m— th order term in Eq. (226) can be evaluated as follows: 

y/hka (4P'^^^VV6n,) = ne^ (4P^*^'^V;,n,) = ne^^ (-4P^o^°F^oo) 

= (4mg^°) . (228) 

where Q'^^^'^ = (l/m)P°'"^'^. Rest of the calculation proceeds exactly as from Eq. (105) 
to Eq. (108) and we find that Eq. (226) gives the horizon entropy. This is a clear reason 
why we can think of S as entropy. 

7.5. Cosmological constant and gravity 

The approach outlined above has important implications for the cosmological constant 
problem [142] which we shall now briefly mention. In the conventional approach, we 
start with an action principle which depends on matter degrees of freedom and the 
metric and vary (i) the matter degrees of freedom to obtain the equations of motion for 
matter and (ii) the metric g"-'^ to obtain the field equations of gravity. The equations 
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of motion for matter remain invariant if one adds a constant, say, — po to the matter 
Lagrangian. However, gravity breaks this symmetry which the matter sector has and po 
appears as a cosmological constant term in the field equations of gravity. If we interpret 
the evidence for dark energy in the universe (see ref. [143]; for a critical look at data, see 
ref. [144] and references therein) as due to the cosmological constant, then its value has 
to be fine-tuned to satisfy the observational constraints. It is not clear why a particular 
parameter in the low energy sector has to be fine-tuned in such a manner. 

In the alternative perspective described here, the functional in Eq. (205) is clearly 
invariant under the shift — Po or equivalently. Tab Tab + PoQab, since 

it only introduces a term —pQUaU"" = for any null vector Ua- In other words, one 
cannot introduce the cosmological constant as a low energy parameter in the action in 
this approach. We saw, however, that the cosmological constant can re-appear as an 
integration constant when the equations are solved. The integration constants which 
appear in a particular solution have a completely different conceptual status compared 
to the parameters which appear in the action describing the theory. It is much less 
troublesome to choose a fine-tuned value for a particular integration constant in the 
theory if observations require us to do so. From this point of view, the cosmological 
constant problem is considerably less severe when we view gravity from the alternative 
perspective. 

This extra symmetry under the shift Tab — ^ Tab + Podab arises because we are not 
treating metric as a dynamical variable in an action principle."'" In fact one can state 
a stronger result [145, 146]. Consider any model of gravity satisfying the following 
three conditions: (1) The metric is varied in a local action to obtain the equations of 
motion. (2) We demand full general covariance of the equations of motion. (3) The 
equations of motion for matter sector are invariant under the addition of a constant to 
the matter Lagrangian. Then, we can prove a 'no-go' theorem that the cosmological 
constant problem cannot be solved in such model. That is, we cannot solve cosmological 
constant problem unless we drop one of these three demands. Of these, we do not want 
to sacrifice general covariance encoded in (2); neither do we have a handle on low energy 
matter Lagrangian so we cannot avoid (3). So the only hope we have is to introduce an 
approach in which gravitational field equations are obtained by varying some degrees 
of freedom other than gab in a maximization principle. This suggests that the so called 
cosmological constant problem has its roots in our misunderstanding of the nature of 
gravity. 

Our approach is not yet developed far enough to predict the value of the 
cosmological constant. But providing a mechanism in which the bulk cosmological 
constant decouples from gravity is a major step forward. It was always thought that 

+ It is sometimes claimed that a spin-2 graviton in the hnear Hmit has to couple to Tat in a universal 
manner, in which case, one will have the graviton coupling to the cosmological constant. In our 
approach, the linearized field equations for the spin-2 graviton field hab = gab — Vat, in a suitable 
gauge, will be {Ohab — Tab)n"Ti'' ~ for all null vectors n". This equation is still invariant under 
Tab ^ Tab + Po9ab showing that the graviton does not couple to cosmological constant. 
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some unknown symmetry should make the cosmological constant (almost) vanish and 
weak (quantum gravitational) effects which break this symmetry could lead to its small 
value. Our approach provides a model which has such symmetry. The small value 
of the observed cosmological constant has to arise from non-perturbative quantum 
gravitational effects at the next order, for which we do not yet have a fully satisfactory 
model. (See, however, Ref. [147, 148].) 

7. 6. Thermodynamic route to gravity: Summary of the paradigm 

The paradigm described in this review is summarized in Fig. 5. The key idea is that the 
behaviour of bulk spacetime is similar to the behaviour of a macroscopic body of, say, gas 
and can be usefully described through thermodynamic concepts — even though these 
concepts may not have any meaning in terms of the true microscopic degrees of freedom. 
This is exactly similar to the fact that, while one cannot attribute entropy, pressure or 
temperature to a single molecule of gas, they are useful quantities to describe the bulk 
behaviour of large number of gas molecules. In such a paradigm, the field equations can 
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be obtained by extremising the entropy expressed in terms of suitable variables. The 
motivation for such a thermodynamic route to gravity is amply demonstrated by the 
existence of local Rindler observers who perceive local horizons and thermal behaviour 
around any event in spacetime and is summarised in the boxes in Fig. 5 leading to the 
central theme: "Spacetime has an entropy density" from the top. 

The form of the entropy function encodes the information about dynamics and its 
extremisation leads to the field equations. We have described the specific forms of this 
function in different contexts along with their physical meaning and inter-relationship 
in the previous sections (see e.g., Eq. (204), Eq. (219) and Eq. (197)). In the context of 
quadratic functionals, we are led to Lanczos-Lovelock model of gravity in general and 
to Einstein gravity (uniquely) in D = 4. This is indicated in the boxes at the bottom 
of Fig. 5. 

In the complete description (statistical mechanics of the 'atoms of spacetime') one 
should be able to obtain the form of the entropy in terms of the microscopic degrees of 
freedom (as indicated by dashed arrows in Fig. 5). In the absence of such a theory, we are 
relying on a thermodynamic description in which the form of the entropy function that 
leads [131] to equipartition of energy among the degrees of freedom and to acceptable 
field equations. The leading term for entropy function in the correct theory is likely to 
be quadratic, thereby giving rise to Lanczos-Lovelock models, but in a full description 
we will also be able to calculate further corrections. 

8. Conclusions and outlook 

This review concentrated, true to the title, several new insights which have been gained 
regarding the thermodynamical aspects of gravity. It presented a case arguing that: (a) 
The conventional approach in which one considers thermodynamical aspects of gravity as 
an interesting but subsidiary results of, say, doing quantum field theory in a spacetime 
with horizon, is fundamentally flawed, (b) There are several features of the theory 
which should be considered as strong hints favouring a fundamental revision of our 
approach towards gravity and spacetime. These hints appear in the form of peculiar 
relationships, especially in the structure of the action functionals describing gravity, 
and in the generality of the thermal phenomena they represent, (c) The last part of the 
review presented an alternate perspective which holds the promise for providing a more 
natural backdrop for understanding the relationship between gravitational dynamics 
and horizon thermodynamics. 

It is useful to distinguish clearly (i) the mathematical results which can be rigorously 
proved from (ii) interpretational ideas which might evolve when our understanding of 
these issues deepen. 

(i) From a purely algebraic point of view, without bringing in any physical 
interpretation or motivation, we can prove the following mathematical results: 

• Consider a functional of null vector fields n°'{x) in an arbitrary spacetime given by 
Eq. (205) [or, more generally, by Eq. (216)]. Demanding that this functional is an 
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extremum for all null vectors leads to the field equations for the background 
geometry given by {2Qab — Tab)n""nf' = where Qab is given by Eq. (97) [or, more 
generally, by Eq. (81)]. Thus field equations in a wide class of theories of gravity can 
be obtained from an extremum principle without varying the metric as a dynamical 
variable. 

• These field equations are invariant under the transformation Tab — > Tab + Podab, 
which relates to the freedom of introducing a cosmological constant as an integration 
constant in the theory. Further, this symmetry forbids the inclusion of a 
cosmological constant term in the variational principle by hand as a low energy 
parameter. That is, we have found a symmetry which makes the bulk cosmological 
constant decouple from the gravity. When linearized around fiat spacetime, the 
graviton inherits this symmetry and does not couple to the cosmological constant. 

• On-shell, the functional in Eq. (205) [or, more generally, by Eq. (216)] contributes 
only on the boundary of the region. When the boundary is a horizon, this terms 
gives precisely the Wald entropy of the theory. 

It is remarkable that one can derive not only Einstein's theory uniquely in D = 4 but 
even Lanczos-Lovelock theory in Z) > 4 from an extremum principle involving the null 
normals without varying gab in an action functional!. 

(ii) To provide a physical picture behind these mathematical results, it is necessary 
to invoke certain effective degrees of freedom which can participate in the (observer 
dependent) thermodynamic interactions near any local patch of a null surface that acts 
as a horizon for certain class of observers. At present we have no deep understanding 
of how this comes about but, at a qualitative level, the physical picture is made of the 
following ingredients: 

• Assume that the spacetime is endowed with certain microscopic degrees of freedom 
capable of exhibiting thermal phenomena. This is just the Boltzmann paradigm: 
// one can heat it, it must have micro structure! \ and one can heat up a spacetime. 

• Whenever a class of observers perceive a horizon, they are "heating up the 
spacetime" and the degrees of freedom close to a horizon participate in a very 
observer dependent thermodynamics. Matter which flows close to the horizon (say, 
within a few Planck lengths of the horizon) transfers energy to these microscopic, 
near-horizon, degrees of freedom as far as the observer who sees the horizon is 
concerned. Just as entropy of a normal system at temperature T will change by 
5E/T when we transfer to it an energy 5E., here also an entropy change will occur. 
(A freely falling observer in the same neighbourhood, of course, will deny all these!) 

• We proved that when the field equations of gravity hold, one can interpret this 
entropy change in a purely geometrical manner involving the Noether current. From 
this point of view, the normals n°- to local patches of null surfaces are related to 
the (unknown) degrees of freedom that can participate in the thermal phenomena 
involving the horizon. These degrees of freedom seem to obey standard rules of 
thermodynamics, including equipartition. 
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• Just as demanding the validity of special relativistic laws with respect to all freely 
falling observers leads to the kinematics of gravity, demanding the local entropy 
balance in terms of the thermodynamic variables, as perceived by local Rindler 
observers, leads to the field equations of gravity in the form {2Qab — Tab)n°'n^ = 0. 

As stressed in earlier sections, this involves a new layer of observer dependent 
thermodynamics. In particular, since observers in different states of motion will have 
different regions of spacetime accessible to them — for example, an observer falling 
into a black hole will not perceive a horizon in the same manner as an observer who is 
orbiting around it — we are forced to accept that the notion of entropy is an observer 
dependent concept. Fundamentally, this is no different from the fact different freely 
falling observers will measure physical quantities differently comapared to non-geodesic 
observers; but in this case standard rules of special relativity allow us to translate the 
results between the observers. We do not yet have a similar set of rules for quantum field 
theory in noninertial frames. It seems necessary to integrate the entire thermodynamic 
machinery (involving what we usually consider to be the 'real' temperature) with this 
notion of LRFs having their own (observer dependent) temperature. 

This requires an intriguing relationship between quantum fluctuations and thermal 
fluctuations. As an illustration, consider the relation 6E = T6S obtained in Sec. 4.4 
for the one-particle excited state which has a curious consequence when we take the 
nonrelativistic limit [131]. The mode function e*"^'^^*'''*(O|0(x)|lk) corresponding to a 
one-particle state in either inertial frame or Rindler frame goes over to a wave function 
i/j{x) in the nonrelativistic (c — >■ oo), quantum mechanical, limit such that ip{x) satisfies 
a Schrodinger equation with an accelerating potential V = mgx when viewed from the 
second frame. In describing the motion of a wave packet corresponding to such a particle, 
the quantum mechanical averages will satisfy the relation (SE) = mg{Sx) = F{6x). On 
the other hand, given the thermal description in the local Rindler frame, we would expect 
a relation like {6E) = TAS to hold suggesting that the entropy gradient AS (due to the 
gradient An in the microscopic degrees of freedom) present over a region {6x) to give rise 
to a force F = TAS/ {6x). If one assumes that (i) AS/ks has to be quantized (based on 
the results of ref.[135]) in units of 2-7r and (ii) {6x) ~ h/mc for a particle of mass m, then 
we reproduce F = mg on using the Rindler temperature ksT = hg/27rc. Alternatively, if 
one assumes that the force F = TAS/ {6x) should be equal to mg, then the universality 
of the Rindler temperature for bodies with different m arises if we use {5x) = h/mc. 
In this case — which involves the quantum mechanical limit of a one-particle state in 
a non-inertial frame — we need to handle simultaneously both quantum and thermal 
fluctuations. The expression F = TAS/ {6x) demands an intriguing interplay between 
thermal fluctuations (in the numerator, TAS, arising from the non-zero temperature 
and entropy in local Rindler frame) and the quantum fluctuations (in the denominator, 
(Sx), related to the intrinsic position uncertainty h/mc) for the theory to be consistent, 
including the choice of numerical factors. 

At a conceptual level, this may be welcome when we note that every key progress in 
physics involved realizing that something we thought as absolute is not absolute. With 
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special relativity it was the flow of time and with general relativity it was the concept 
of global inertial frames and when we brought in quantum fields in curved spacetime it 
was the notion of particles and temperature. 

Many of these technical issues possibly can be tackled in more or less straightforward 
manner, though the mathematics can be fairly involved. But they may not be crucial 
to the alternative perspective or its further progress. The latter will depend on more 
serious conceptual issues, some of which are the following: 

(i) How come the microstructure of spacetime exhibits itself indirectly through 
the horizon temperature even at scales much larger than Planck length and obeys an 
equipartition law (see e.g. Eq. (180) and ref. [69, 131])? This is possibly because 
the event horizon works as some kind of magnifying glass allowing us to probe trans- 
Planckian physics [82, 83] but this notion needs to be made more precise. 

(ii) How does one obtain the expression for spacetime entropy density from some 
microscopic model? In particular, such an analysis — even with a toy model — should 
throw more light on why normals to local patches of null surfaces play such a crucial role 
as effective degrees of freedom in the long wavelength limit. Of course, such a model 
should also determine the expression for P"-^'^'^ and get the metric tensor and spacetime 
as derived concepts - a fairly tall order!. (This is somewhat like obtaining theory of 
elasticity starting from a microscopic model for a solid, which, incidentally, is not a 
simple task either.) 
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